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ABSTRACT 


This  finsl  rsport  describes  the  work  perforaed  on 
various  coabustion  probieas  related  to  high  frequency 
instability  in  liquid  rocket  engines.  Using  the 
steady-stat*  and  instability  coaputer  prograas  developed 
under  this  study,  a  paraaetric  investigation  was  con¬ 
ducted.  This  investigation  deterained  the  influence  of 
droplet  radius,  droplet  distribution,  injection  velocity, 
chaaber  pressure,  and  aixture  ratio  on  the  ainiaua  threshold 
disturbance  required  to  trigger  coabustion  instability  in 
a  Transtage  type  engine  configuration.  The  propellant  con- 
bination  considered  was  nonoaethyihydrazine/nitrogen  tetro- 
xide.  Results  of  the  study  show  that  increases  in  injection 
velocity  and  droplet  distribution  increased  stability. 

An  increase  in  chaaber  pressure,  based  on  constant  flow 
rate,  increased  stability  while  increases  in  chaaber 
pressure,  at  a  constant  contraction  ratio,  decreased 
stability.  There  appears  to  be  a  droplet  size  for 
ainiaua  stability,  with  changes  in  either  direction 
resulting  in  iaproved  stability.  Results  also  show  that 
due  to  the  vapor  phase  reactions,  aonoaethylhydrazine/ 
nitrogen  tetroxide  vaporize  at  approxiaately  the  sane 
rate.  Thus,  the  oxidizer  or  fuel  could  be  aade  to 
control  the  coabustion  process  by  slight  changes  in  the 
injector  and  engine  paraaeters.  For  the  engine  configura¬ 
tion  studied  the  oxidizer  vaporized  slower  thsn  the  fuel. 

A  prograa  review,  work  on  droplet  atoaization,  and  a 
study  of  hypergolic  liquid  phase  reactions  is  also 
reported.  The  review  suaaarizes  the  work  reported  in 
the  Seaiannual  (DSC  SN-68-51)  and  Special  report 
(AFRPL-TR-6S-254)  written  under  this  prograa. 
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I.  INTRODUCTION 


The  work  conducted  under  AF  04(61 l)-iOS42  hee  been  directed 
toward  developing  nethoda  for  deaigning  stable  rocket  engines. 
The  various  areas  investigated  during  this  progran  are: 

1)  A  steady-state  conbustion  nodel  for  hydrasine  type 
propellants. 

2)  A  conbustion  instability  nodel. 

3)  A  scaling  criteria  for  stable  engines. 

4)  An  experinental  progran  to  verify  the  developed 
scaling  criteria. 

5)  A  survey  of  atonisation  literature. 

6)  An  experinental  progran  ta  determine  the 
influence  of  propellant  properties  and  hydraulic 
paraneters  on  atonisation. 

7)  An  analysis  of  the  effects  of  solid  particles 
on  danping  conbustion  instability  in  gelled 
propellants. 

8)  A  parametric  study  of  injector  and  engine 
paraneters  on  stability  of  a  Transtage  type 
engine. 

9)  A  study  to  determine  the  effects  of  liquid 
phase  reactions  on  atomisation  and  combustion 
of  hypergolic  propellants. 

10)  A  correlation  and  prediction  of  droplet  site 
from  an  impact  energy  method  by  minimizing 
surface  energy. 

Items  (1)  and  (2),  i.e.,  the  development  of  the  steady-state 
and  combustion  instability  models,  their  application  to  pre¬ 
dict  instability  zones  and  listings  of  the  developed  computer 
programs  are  repo~*ted  in  the  Special  Report  (Ref.l)  and 
reviewed  in  Section  II.,  Items  (3)  to  (7)  were  reported 
in  the  Semiannual  Report  (Ref. 2)  and  pertinent  results  are 
also  summarized  in  Section  II.  In  Sections  HI,  IV,  and  V, 
Items  (8),  (9),  and  (10)  will  be  covered  respectively.  Thus, 
this  report  vill  concentrate  on  the  results  of  a  parametric 
study  on  conbustion  instability  ir  a  Transtage  type  engine, 
the  effects  of  liquid  phase  reactions  on  atomization  and 
combustion,  and  a  method  of  predicting  droplet  size  from 
an  impact  energy  method  by  minimizing  surface  energy. 
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II.  REVIEW 


This  ssction  will  revitw  the  work  conducted  on  the 
subject  contract  and  reported  in  detail  in  the  Special 
and  Seaiannual  Report  (Refs.  1  and  2).  This  section  will 
discuss:  1)  developaent  of  the  steady-state  nodel;  2)  devel- 
opaent  of  a  coabustion  instability  aodel;  3)  criteria  for 
scaling  and  an  erp*'.iaental  program  to  verify  these  criteria 
4)  an  atoaisation  literature  review  and  an  outline  of  an 
experimental  program  to  determine  the  influences  of  pro¬ 
pellant  parameters  and  triplet  injector  geometries  on  drop¬ 
let  size;  and  S)  an  analysis  of  the  effects  of  solid  particles 
on  damping  combustion  instability  in  gelled  propellants. 

Steady-State  Combustion  Model/Hydrazine-WTO 

Analysis  of  high  frequency  combustion  instability 
requires  complete  spatial  knowledge  of  such  parameters 
as  vaporization  rate,  relative  velocity,  and  droplet 
Reynolds  number.  To  determine  these  parameters,  a 
detailed  steady-state  combustion  model  was  developed 
for  decomposing  type  propellants  and  additAonal  informa¬ 
tion  was  determined  on  hydrazine  droplet  combustion, 
atomization  and  liquid  phase  reactions.  To  develop  a 
hydrazine  combustion  model,  experiments  were  conducted 
which  indicated  that  a  single  diffusion  flame  front 
model,  as  shown  in  Figure  1,  is  not  realistic  and  two- 
flame  fronts  appear  in  the  combustion  of  MMU  and  hy¬ 
drazine  in  nitrogen  tetroxide,  as  shown  in  Figures  2 
and  3.  It  has  been  shown  that  existence  of  "two  flames" 
depends  on  the  convective  environment  and  the  hydrazine 
decomposition  kinetics.  The  two-regime  model  of  the 
decomposition  front  position  is  shown  in  Figure  4.  In 
regime  I  the  laminar  diffusion  mantle  thickness  B  is  • 
larger  than  the  distance  required  for  decomposition( X*) 
and  two  flames  occur  simultaneously,  i.e. ,  an  outer 
oxidation  flame  with  an  inner  decomposition  flame. 

Regime  II  is  the  same  as  the  laminar  diffusion  model 
shown  in  Figure  1.  Here  the  film  thickness  is  thinner 
than  the  distance  required  for  the  hydrazine  to  decompose 
and  a  single  flame  exists.  Profiles  of  combustion  rate 
and  film  thickness  along  the  combustion  chamber  length 
are  shown  in  Figure  5.  When  the  propellant  is  injected 
into  the  combustion  chamber  it  is  traveling  at  approxi¬ 
mately  the  injection  velocity,  while  the  axial  velocity 
of  the  combustion  gases  is  approximately  zero.  As  com¬ 
bustion  takes  place  the  gas  velocity  increases,  while 


2 


determined  from  Nusselt  number  correlation 


FIGURE  1.  LAMINAR  DIFFUSION  MODEL,  SINGLE  FLAME  FRONT 


MMH  Burning  in  still  nitrogen 
tetroxide,  at  one  atmosphere 


FIGURE  2. 


PHOTOGRAPH  OF  EXPERIMENTAL  BURNING 
DROP  SHOWING  TWO  FLAME  FRONTS 
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Note;  The  liquid 
drop  is  within  the 
wire  support. 


Hydrazine  burning  in  still  nitrogen 
tetroxide,  at  one  atmosphere 


FIGURE  3.  PHOTOGRAPH  OF  EXPERIMENTAL 

BURNING  DROP  SHOWING  TWO  FLAME 
FRONTS 


5 


vAi»r9  B  depends  upon  NiiseeU  nuiriber  of 
correlation 


DBC01IP081TION  FRONT  WITHIN 
DIFFUSION  MANTLE,  REGIME  I 


dtifuelon  flame  front 
together  when  a*  ^  B 
Note: 

X*  is  independent  of 
flow  while  B  decreases 
with  velocity 


\ 


laminar  diffusion  model  when 

X*  ^B,  REGIME  n  (This  is  the  same 
as  Figure  1 ) 


FIGURE  4.  TWO'REGIME  MODEL  OF  DECOMPOSITION 
FRONT  P08ITS0N 
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Film  thickness*  B  (microns)  Combustion  rate*  reCfraction/inclO 


0. 


Chamber  length,  inches 


FIGURE  5.  Combustion  Rate,  Initial  Drop  Size  49.  4y 
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t%m  ti^«i4  intil  tli«  gst  wm4  arc 

travaliag  at  tM  aaaa  valaeity*  Aa  caaAaitiaa  graaaa4t« 
the  faa  ralaaitjr  iacraaaaa  fartkart  vltli  the  liq«i4 
accalaratiag  tairar4  the  gaa  vaiacity*  Tha  eaavaativa 
fils  thiahaaaa  variaa  iavaraaiy  ta  tha  ralativa  valaeity 
hatwaaa  tha  gaaaa  aa4  liatti4«  Praa  Figwre  S  it  is  shawa 
that  whaa  tha  prapallaat  is  first  iajaata4  iato  tha  chaabar 
tha  fila  thichaass  is  lass  thaa  tha  dacaapasitiea  thick* 
aass  hacaasa  af  tha  high  ralativa  valaeity  aad  tha  4raplcts 
ara  haraiag  ia  ragiaa  li«  Tha  gas*liqui4  ralativa  valaeity 
4aeraasas  ta  paiat  k,  ahara  tha  axidatiaa  aa4  Aacaapasitiaa 
ara  tha  saaa«  Up  ta  this  paiat  tha  caahastiaa  rata  fa Haas 
tha  siagla  flaaa  aa4al  salutiaa  as  shaaa  at  paiat  B.  As  the 
ralativa  valaeity  is  further  4acraasa4*  tha  file  thickaass 
graus  larger  thaa  tha  daeaapasitiaa  thickaass  aad  the 
ragiaa  1  eaatrals  tha  eaidiustiaa  praeass.  Tha  iacrassad 
rata  af  tha  **tva-flsaa**  aadal  aecurs  at  paiat  C.  After 
tha  ralativa  valaeity  reaches  a  aialauBf  tha  gas  passes 
tha  draplats  aad  the  fila  thickaass  daeraasas  ta  paiat  E» 
uhara  tha  draplats  ra-aatar  ragiaa  1*  Tha  huraiag  rata 
raaaias  relatively  eaastaat  batwaaa  C  aad  D  siaca  tha 
draplats  have  heated  ta  their  vat^bulh  paiat.  Tha  lata* 
grated  result  af  tha  iacraasad  caabustiaa  rata  is  that 
tha  fraetiaa  af  prapallaat  eaabustad  aaar  tha  iajactar 
faea  is  higher  thaa  uauld  accur  with  prapallaat  uadar 
tha  saaa  eaaditiaas«  without  dacoaposition. 

Ta  eaaputa  staady*stata  caabustiaa  prefixes  a 
detailed  kaawladga  af  tha  vaparixatiaa  rata  of  N2O4  is 
alsa  required.  As  shawa  ia  Figure  tha  aatariag 
axiditar  draplat  is  assuaad  to  be  N2O4.  However,  as  tha 
N2O4  evaporates  it  goes  to  approxiaataly  S0%  NO2.  Than, 
siaca  tha  taaparatura  gradient  is  vary  steep,  (Figure.  6) 
tha  vapor  will  be  100%  NO,  before  it  crosses  even  10%  of 
the  diffusion  aaatla.  The  approxiaation  that  NO2  is  tha 
vaporising  species  was  thus  aade.  Using  this  approxiaation 
it  was  possible  to  account  for  chenical  change  by  adding 
the  Jiaat  of  dacoaposition  to  the  heat  of  vaporisation 
for  N2O4.  Since  this  reaction  is  endotheraic  the  effec¬ 
tive  heat  of  vaporisation  is  increased  and  the  vaporisa¬ 
tion  rate  is  slower  thaa  that  predicted  when  dacoaposition 
is  not  taken  iato  account. 

A  propellant  vaporisation  aodel  developed  by 
Priea  aad  Haidnann  (Ref.  3)  was  used  to  coapare  calcu¬ 
lated  aad  experiaental  coabustion  efficiencies  for  a 
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TIm  tMamptioa  poMibUi  tte  lucliisloii  oC 

beat  of  reacUon  with  heat  of  TaportaatiOB  la 

At  <10%  of  (ft  -r  ). 


DROPLIT  VAraUZATION  MODEL 


FIGURE  6. 
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imrgt  maitt  0i  ftfllmt  CMlbiaatiMis*  (ik«im  itf  FigiiM 
T|«  Miil«  •xcailMt  ««•  with  May 

fyaf«llMts«  afficiMciat  ««r«  alwayi  calealated 

tawar  tbaa  aaaaafa4«  Haiag  tkm  hjrJraa^a  *'t«a-flaaa** 
4acaafaaitiaa  aa4al«  caofcaatiaa  afficiaaciaa  far  iiy4r«- 
xiaa  aa4  gasaaM  asygaa  aara  raeai^atag  as  aliawa  ia 
figara  7«  Tha  raaalta  sha«a4  that  tha  ^acaagaaitian 
■a4ai  fra4icta4  a  aara  accarata  aa4  highar  caahactioa 
rata  thaa  tha  aiagla  fiaaa  aagal  aaa4  hy  Friaa  aa4 
Maidaaaa  (iaf«  1) «earralatiag  tha  4ata«  It  «aa  also 
ahaaa  that  tha  aaaiaaa  affact  acc«rra4  far  shartar 
chaahar  laagtha«  agraaiag  with  tha  4acaii9aaitiaa  ao4ai 
trhich  ahaaa4  tha  iacraaaa  ia  caahaatiaa  aaar  tha  iajac- 
tar  faca  at  tha  aiaiaaa  raiativa  valaaity  ragiaa. 

Ta  ahtaia  a  caaglata  spatial  kaatrlagga  of 
r^ariaatiaa  rata«  raiativa  valacity  aa4  4roplat  Rayael4a 
aiMriiar  tha  ataa4y*atata  vaparixatiaa  aa4al  vith  a  vapor 
phaaa  gacaa^oaitioa  fiaaa  oaa  prograaaag  aa4  aolva4  oa 
a  Caatral  Oata  3444  4igital  coi^wtar^  Coapatar  aolatioat 
ara  haiag  «aa4  ia  paraaatric  ata4iaa  to  4ataraiaa  tha 
affact  af  physical  propartias  aa4  4asiga  paraaatara 
apaa  caaiiastioa#  Tha  coapotar  ao4al  tracts  4acoaposition 
hath  ia  tha  aitragaa  tatroxi4a  aa4  hy4rasiaa  systaa. 

Tha  pragraa  caa  coapiata  vaporixatiaa  of  a  bipropallaat 
systaa«  thas  aaahliag  a  raalistic  pra4ictioa  of  chaahar 
paraaatara. 

2.  lastahility  Mo4al 

Oyaaaic  Sciaaca  Corporatioa  has  4avalopa4  a 
aoaliaaar  ao4al  for  dataraiaiag  tha  xoaas  of  a  liquid 
rockat  aagiaa  ia  vhich  a  taagaatial  aodo  of  high  fra- 
qaaaey  iastahility  is  aost  aasily  iaitiatad.  This  aodal 
has  baaa  ralatad  to  hardvara  dasiga  paraaatars.  i.a., 
paraaatara  ralatad  to  iajactor  dasiga,  chaabar  coa- 
figuratioa,  aad  propallaats,  tharaby  aaabliag  the 
iafluaacas  of  s^xtaa  dssiga  aad  stability  ratiag  davices 
to  ba  dataraiaad. 

Coabustioa  iastiibjklity  ?9nas  ia  a  liquid  rockat 
aagiaa  ara  dataraiaad  by  coabioiag  &  aoaHaaar  iasta- 
bility  aodal  with  a  staady-stata  vaporizs:ioa  prograa* 

Tha  iastability  aodal  coasidars  the  aoaliaaar  coasarva- 
tioa  aquatioas  with  aass  additioa  usiag  a  staady-stata 
axprassioa  for  tha  buraing  rata.  This  aodal  applies  to 
a  ona-diaansioual  annulus  of  saall  length  (iz)  and 


10 


thickntfs  (Ar)  ihown  in  Figurt  g.  Applying  tht 
retuitt  of  this  aodol,  ns  shown  in  Figuro  9,  s  rocket 
ongino  cnn  be  nnalysed  by  increaentnily  dividing  the. 
coabustion  chsaber  into  nnnulnr  nodes  in  the  r  end  t 
directions.  Stendy-stnte  properties  nt  etch  tnnuinr 
node  or  position  in  the  chtaber  ere  computed  from  t 
stendy-stnte  vnporizntion  program  described  in  Fnrngrnjph 
1.  These  steady-state  properties*  i.e**  vaporization 
rate,  relative  velocity,  and  droplet  Reynolds  number, 
and  the  curves  from  the  instability  model  (Figure  9) 
are  used  to  determine  the  stability  of  the  node*  This 
process  is  repeated  for  each  node  to  determine  a  sta¬ 
bility  map  of  the  entire  engine,  as  shown  in  Figure  10, 

From  a  nonlinear  model  important  nonlinear 
phenomena  are  predicted,  i.e.,  (I)  stability  dependence 
on  disturbance  wave  shape,  amplitude,  type  (velocity 
or  pressure),  and  position;  (2)  the  limiting  amplitude 
of  the  unstable  pressure  oscillations;  and  (3)  the  shape 
of  the  unstable  wave  forms.  Major  results  of  this  analy&i 
show  that  the  amplitude  and  position  of  a  pressure  dis¬ 
turbance  required  to  initiate  instability  can  be  deter¬ 
mined,  thereby  defining  a  sensitive  zone  (and  the  best 
plaice  to  disturb  vhe  engine)*  This  sensitive  zone 
extends  severs  .,  f.r  ches  from  the  injector  face  and  occurs 
where  the  tiyr&gj  droplets  are  moving  the  slowest,  relative 
to  the  gases* 

Dynamic  Science  Corpuration  has  shown  that  an 
annular  combustor  section  will  be  more  stable  as  the 
droplet  Reynolds  number  approaches  zero*  Thi.^s,  a  signi¬ 
ficant  result  of  this  work  is  that  there  are,  i4x  a 
vaporization  controlled  combustion  process,  three 
parameters  affecting  stability: 

(a)  Buring  rate  parameter  -  L 

(b)  Absolute  value  of  relative  velocity  -  Av 

(c)  Reynolds  number  of  drop  based  on  speed 
of  sound  •  Re^. 

Such  information  is  not  only  of  use  from  the 
preliminary  design  standpoint,  but  also  as  an  invaluable 
tool  in  understanding  how  rocket  engines  should  be 
disturbed  during  their  development  test  programs  to 
determine  the  degree  of  stability  of  an  engine.  The 
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FIGURE  10.  Pulse  Pressure  Contour  Lines  Defining  Unstable  Zones 
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Ofstaic  Scitact  Corporatioa  aodal  taablet  tha  tngiae 
datigaar  to  dlataraiaa  tka  poaitioa  to  iatrodoca  tha 
dittarbaaca,  a  raasonabla  dittarbaaca  aaplitada 
critaria,  aad  tha  aoat  affactiva  wava  prof i la  (prastura 
aad/or  valocity  diiturbaaca) .  Siaca  tha  iajactor  dasign 
variablaa  caa  ba  ralatad  to  thrashold  disturbaaca 
aaplituda,  paraaatara  caa  ba  aodifiad  to  iacraasa 
stability  of  a  givaa  aagiaa  coafigaration  or  ba 
ttsad  ia  tha  praliaiaary  dasiga  of  aa  aagiaa. 
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Scaling. 

Otsign  anU  scaling  criteria  for  stable  liquid  rocket  engines 
were  developed  using  a  nonlinear  coabustion  instability  aodel  with 
a  droplet  vaporisation  aodel.  The  scaling  paraaeters  were  foraed 
into  an  ''Instability  Nuaber"  which  was  used  to  correlate  approxi- 
aately  60  runs  conducted  by  the  USAF  RPL  and  Aerojet. 

Result  of  the  nonlinear  instability  showed  that  three  para- 
■aters  define  the  shape  of  a  sensitive  instability  zone:  (1)  burning 
rate  paraaeter;  (2)  Reynolds  nuaber  of  the  droplets  based  on  the 
speed  of  sound;  and  (3)  relative  velocity  between  the  gas  and 
droplets.  The  study  further  showed  that  a  liquid  rocket  engine 
is  nost  sensitive  when  the  gas-droplet  relative  velocity  is  a 
ainiaun  or  zero.  Thus«  in  developing  a  scaling  paraaeter  for 
the  tangential  node  of  high  frequency*  one  of  the  three  inportant 
paraneters  can  be  eliainated  and  only  the  burning  rate  paraaeter 
(L)  and  (Re^)  Reynolds  nuaber  of  droplet*  need  be  considered  at 
6V«0.  If  two  engines  have  the  sane  L  and  Re^  at  AV»0*  they  will 
have  the  sane  stability*  that  is*  the  sane  ainiaun  threshold 
disturbance  required  to  trigger  instability. 

For  nost  injectors  built*  the  aost  sensitive  region  is 
at  the  outer  radius  of  the  engine.  Therefore*  if  there  is  a 
snail  dianeter  low- thrust  high-perforaance  engine  utilizing  a 
given  propellant*  and  a  larger  thrust  engine  with  the  sane 
stability  characteristics  is  required*  scaling  the  burning  rate 
paraaeter  and  Reynolds  nuaber  of  the  droplet  at  the  outer  radius 
on  the  AV«0  plane  will  insure  the  sane  stability  characteristics 
(for  L<2).  However*  since  the  injector  variables  will  be  changed 
in  this  scaling  process*  length  as  well  as  dianeter  aust  be 
scaled  to  retain  n^«.  As  shown  in  Figure  11*  scaling  froa  a  nodal 
engine  (1)  to  a  full-size  engine  (2)  required  (lAV-O  and  L<2) : . 

•^•di  i  ‘^®d2 

Li  >  L2 

By  scaling  with  conditions  (11-1)  and  (17-2)  the  full-size  engine 
has  equal  or  greater  stability*  i.e.*  an  equal  or  greater  disturbance 
required  to  trigger  instability.  By  assuming  that  the  propellant 
and  aii^ture  ratio  are  retained  in  the  scale  engine  L  and  Re^  nay 
be  related  to  engine  paraneters  by 

.  iLh  <  1 

"•d,  rj  Pj  - 
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Disturbance  Pressure  Amplitude,  AP/^ 


h 
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>^ini  ?!  ^2 
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Th«s«  criteria  are  ahoMn  in  Figure  12.  For  conditions  11*1 
and  11-2  to  be  fulfilled  the  ratios  aust  fall  inside  the  box*  To 
coabine  these  two  groups  a  functional  relationship  is  required 
between  L,  Re^  and  AP/P  such  that: 


For  L<2,  such  a  relationship  can  be  approxiaated  by  curve  fitting 
the  coaputed  stability  liait  curves,  with  an  expression  of  the  fora: 

AP  1 

a  , . . -  where  n«l  to  3 

^  L  Re^n 

Lines  1  and  2  of  Figure  12  represent  locus  of  solutions  for  n«l 
and  3  respectively.  The  conditions  Il-l  and  IX-2  guarantees  (L<2): 


a. 


The  engine  is  aore  stable  than  the  scale  if 


*^*d2 


<1 


b.  The  engine  is  as  stable  as  the  scale  if 

.  S  ;  **<»2  .  I 

•■1  TS37 

c.  The  engine  is  less  stable  than  the  scale  if 


Red2 


>  1 


To  facilitate  experiaental  correlation,  L  and  Re^  are 
foraed  into  an  "instability  nuaber."  Engines  with  the  saae 
instability  nuaber  will  have  the  saae  stability,  i.e.,  they 
will  be  triggered  unstable  to  the  saae  AP/P  pulse.  As  presented 
in  Reference  2^two  instability  nuabers  are  foraed  froa  the 
product  of  Reynolds  nuaber  of  the  droplet  and  burning  rate 
paraaeter. 
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riGURS  .^*  STABILITY  PARAMETERS  (l^<2) 
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Vliere  E  is  the  stuuimtioa  of  the  iapect  energy  which  the  strems 
exert  nornnl  to  each  other.  This  tera  is  an  atteapt  to  relate 
streaa  geoaetry  and  velocity  to  drop  radius. 

These  "instability  nnabers"  are  used  to  correlate  engine 
tests  as  shown  in  Figures  1 3*! 7.  In  addition  to  the  NNH  tests 
conducted  at  USAF-RPL,  this  relation  correlated  S7  out  of  61 
tests  with  pulse  aotors  at  Aerojet*  using: 

a.  MON- 10  and  U-DETA 

b.  N20^  and  N^H^ 

c.  XRFNA  and  UDNH 


Figures  IS-IS  show  the  usual  stability  decreasing  with  increasing 
I  ;  however  at  a  certain  value  of  I  the  trend  reverses  itself 
corresponding  to  aoving  to  the  L>2  side  of  the  stability  curves. 
Pulse  strengths  froa  the  RPL  data  are  lower  because  they  are 
neasured  in  a  rocket  engine,  while  the  Aerojet*  data  was  aeasured 
in  a  shock  tube.  The  Aerojet  data  does  not  coapare  quanitatively 
with  the  RPL  data  due  to  difference  in  aensureaents ,  however  the 
trends  are  the  saae. 


In  conclusion,  a  set  of  scaling  criteria  have  been  developed 
which  enable  retaining  scale  stability  when  building  a  larger  engine 
or  changing  operation  conditions.  An  instability  index  or  nuaber 
is  being  derived  which  enable  rating  engines  in  order  of  Stability. 
Several  seal-analytical  expressions,  based  on  the  product  of  Re. 
and  L  were  successful  in  the  correlation  of  the  experiaental  data 
however,  aore  work  has  to  be  done  before  a  coaplete  understanding 
of  the  stability  nuaber  can  be  obtained. 


* (Sray,  P.D., and  Krieg,  H.C.  "Pulse-Motor  Evaluation  of 
Injector-Pattern  Coabustion  Stability  with  Storable  Propellants, 
Aerojet  Special  Report  for  USAF-RPL,  Confidential,  AO-318667, 

29  July  1960. 


FIGURE  14. 


Instability  Number  Correlation  for  MMH/NTO  Tests 
Pcrforaed  by  U.S.Air  Force  Rocket  Propulsion  Lab 


FIGURE  15. 


Instability  Nuaber  Correlation  for  MON>10  and 
U-DETA  Tests  Perforaed  by  Aerojet. 


FIGURE  16. 


Instability  Nuiiber  Correlation  for  N-0./N„H 
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Tests  Performed  by  Aerojet. 


4,  Droplet  Atoaitation, 


Nuaerous  experiaents  directed  toward  the  specification 
of  size  distribution  of  liquid  droplets  after  atoaization 
have  been  reported  in  the  literature.  Statistical  distribution 
functions  have  been  developed  eapirically  by  fitting  experi¬ 
mental  data.  Spray  distribution  has  been  related  to  both 
physical  properties  and  injector  design  variables;  however, 
it  is  not  possible  to  completely  predict  a  realistic  distri¬ 
bution  and  determine  its  statistical  parameters  based  on  given 
propellant  properties  and  injector  design. 

Actual  drop  distributions  are  not  uniform  across  the  chamber 
but  have  three-dimensional  spatial  characteristics.  The 
velocity  of  droplets  is  also  distributed  over  a  velocity 
range.  The  initial  droplet  velocity  is  as  important  as  the 
initial  size.  However,  information  concerning  droplet  velocity 
distribution  is  current  not  available. 

The  distribution  of  droplets  varies  with  the  mechanisms 
of  atomization.  Two  generalized  categories  nay  be  used  to 
describe  the  wide  range  of  liquid  atomization  element  de¬ 
signs.  These  are  plain-jet  injectors,  and  sheet  atomizers. 

The  most  commonly  used  rocket  injector,  the  impinging  jet, 
combines  the  characteristics  of  the  sheet  and  plain-jet 
injectors.  Ueidmann  and  Foster  (Ref. 4)  and  Bittker  (Ref.S) 
have  suggested  that  the  impinging- jet  injector  produces  a 
binodal  distribution  of  particles,  this  being  consistent  with 
the  hybrid  nature  of  the  impinging- jet  atomization  mechanisms. 

Many  statistical  distribution  functions  have  been  empirically 
fitted  to  experimental  spray  data.  Bevans  (Ref.  6)  has 
examined  the  three  most  commonly  used  distribution  functions, 
the  Nukiyama-Tanasawa,  Rosin-Rammler,  and  logarithmico-normal 
and  concluded  that  the  Rosin-Rammler  expression  was  most 
successful  at  expressing  experimental  data  from  a  pressure- 
atomized  oil  spray.  However,  Ingebo  (Ref.  7,  8,  and  9)  has 
consistently  obtained  good  results  in  fitting  his  excellent 
experimental  work  to  a  Nukiyama-Tanasawa  distribution.  Priem 
and  Heidmann  (Ref.  3)  state  that,  "No  criterion  has  been 
established  for  choosing  the  preferred  distribution  function. 
For  convenience,  therefore,  the  logarithmiconormal  distribution 
function  was  used  in  this  investigation...."  General 
mathematical  forms  of  these  three  distributions  are: 
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Rosin-Rftaaler ; 


iih«r« 

and 


whara 
and 
whi  la 


R  •  axp  (-4/t^)*R 

R  is  tha  aass  fr.tction  of  drops  larger  than  r 

r^i  and  a^  art  constants  which  detarnine  tha  aadian 
site  and  disparsiotk  oi  the  droplets. 

Lcgaxithaico- normal: 

R  is  tha  nass  fraction  of  drops  larger  than  r, 
a|^  is  an  eapirical  constant 

r^  and  Og  are  the  nass  aedian  drop  radius  and  the 
standard  deviation. 


Nuk iyaaa- Tanas awa: 


dR 

dr 


where  a^  is  an  eapirical  constant. 

Tha  values  of  tha  eapirical  constants  shown  in  the  above 
expressions  have  bean  deternined  in  Refs.  4,  7,  8,  and  9 
for  specific  fuels  in  various  injectors.  All  of  these  dis¬ 
tributions  exhibit  large  errors  when  used  in  the  snail  drop 
size  range.  More  experiaental  data  are  required  in  this 
range  so  that  a  preferred  distribution  nay  be  deternined. 


Ingebo  and  Foster  (Ref.  7 )  used  diaensional  analysis  to 
study  the  effect  of  seven  physical  properties  upon  the  volune- 
aean  droplet  radius. 


They  considered  the  following  variables  as  having  a  possible 
influence, 


^30  "fCDo.PU  Oj,  Ui,  p^,  Pg) 
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where:  ■  injector-orifice  dieneter 

u  «  free  streea  velocity 
0  ■  surface  tension 

y  ■  viscosity 
p  ■  density 

After  diaensional  grouping  of  terms  they  fit  their  nuaerous 
data  with  the  equation 


or 
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where: 


-e) 


1/h 


0l 


0  p 

o  ^g 


■  Weber  Nuaber»  We 


and 

Dj,  Pi  u 

-  - -  ■  Reynolds  Nuabe^,  Re 

This  data  was  recorded  for  cross-current  break-up  of  liquid 
jets  in  high  velocity  airstreaas. 


They  also  eapirically  deterained  the  expression  for  the  aaxi 
mua  droplet  radius, 


11.1 


We  \  0.29 

^  / 


and  fit  their  size  population  with  the  Nuki/ama-Tanasawa 
dist  ributi  )n : 


Pri«a  and  lleldaann  (Raf.  3)  used  thast  rasults  (Raf.  7) 
in  tha  fora 


*■30  •  r,  a 

Pi 

tehcra  it  tha  aass  aadian  droplet. 

Howevar,  if  the  standard  deviation  changes,  tha  voluaeoaean  drop 
r.Q*  is  not  proportional  to  r  .  the  aass-aedian  drop.  Priea  and 
Heidaann  used  this  relation  afthough  there  is  no  axperiaental 
proof  of  its  validity.  Using  this  fora  they  coapared  the  aedian 
site  for  various  propellants  with  tha  data  of  Ingabo  (Ref.  8)  for 
n>haptane  iapinging> jet  braak>up  in  airstreaas  siaulating  the 
velocity  conditions  in  rocket  coabustors.  The  coaparison 
equation  was  •  where:  x  is  the  unknown  propellant 

h  is  heptane 


Herrington  and  Richardson  (Ref.  10)  plotted  data  showing  a  varia¬ 
tion  of  four  orders  of  aagnitude  in  y/p .  These  data,  shown  in 
Figure  18.  deteraina  that  r  o(y  /p  )  ;  this  is  very  close  to 

tha  0.25  power  of  Ref.  8  ano  can^be^used  with  confidence.  The 
surface  tension  of  aost  of  the  liquids  of  Figure  18  was  30  dynes/ca 
except  for  water  and  glycerine  which  have  valuer  of  73  and  63 
respectively.  These  two  liquids  produced  aean  drop  sizes 
expected  froa  their  respective  kineaatic  viscosities .  and  showed 
no  surface  tension  effects.  Thus,  the  effect  of  surface  tension 
on  aean  drop  size  is  not  fully  understood. 

Griffen  and  Laab  (Ref.  11)  studied  the  dependence  of  aean  drop- 
size  on  aabient  gas  density.  They  deterained  that  the  aaxiaum 
droplet  size  was  reduced  with  increasing  ambient  gas  density 
according  to  the  relation 


They  coanent  that  the  apparent  iaprovement  in  atomization  (mean 
drop-size)  with  increasing  pressure  is  mainly  due  to  a  reduction 
in  the  number  of  large  drops  and  is  probably  >'aused  by  subdivision 
of  these  big  droplets. 
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Ing«bo  (Ref.  t)  correlated  date  obtained  for  heptane 
sprays  produced  by  pairs  of  90*  inpinging  jets  over  wide 
ranges  of  orifice  diaaeter,  D^;  TTquid-jet  velocity,  v. ; 
and  jet-airstreaa  velocity  difference,  Av;  according  to 
the  equation 

^0  ■  [“•  ’'j]  “o 


Ueidaann  and  Foster  (Ref.  4)  studied  the  effects  of  inpingeaent 
angles  of  10*  to  90*  and  jet  velocities  of  50  to  74  feet  per 
second.  All  their  distributions  showed  biaodal  characteristics. 
The  nost  significant  effect  of  inpingeaent  angle  and  jet  velo* 
city  was  a  change  in  the  relative  nunber  of  drops  in  each 
node.  This  is  consistent  with  the  hybrid  nature  of  inpinging 
jets  (i.e.,  drag  instability  break-up)  and  the  fact  that  they 
deterained  that  sheet  thickness,  t,  varied  according  to 

t  ,  (P 

This  data  (Ref.  4)  could  be  further  analyzed,  using  #n 
instability  nodal  to  relate  sheet  thickness  to  drop  size  dis¬ 
tribution.  Thus  a  single  aodel  and  seni-eapirical  equation 
would  be  developed,  using  injector  paraneters  to  predict 
sheet  thickness,  and  physical  properties  to  predict  insta¬ 
bility  break-up  and  final  size  distribution. 

In  the  absence  of  any  general  theoretical  treataent  of  the 
effect  of  aass  flow  rate  and  velocity  upon  r3o,  several 
purely  eapirical  relations  have  been  deterained  for  different 
atoaiser  designs.  Radcliffe  (Ref.  12)  deterained  the  sur¬ 
face  aean  radius,  T2q,  of  the  spray  produced  by  a  simple 
swirl  atoaizer.  He  developed  the  relation 

rjo  «  /  ApO-530 
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fuel  flow  rate 


>  •:  I  c  h 

and  AP  •  injector  pressure  drop. 

Fraser  and  Eisenklaa  (Ref.  13)  obtained  the  following  correla¬ 
tion  for  a  fan  spray  nozile 

rjo  . 

For  a  soaewhat  different  type  fan  spray  (Ref,  14  and  15)  it 
was  deterained  that 

rjo  .  K0.25/4p0.305 

It  can  be  concluded  that  while  there  are  eapirical  correlations 
for  soae  special  injectors,  it  is  not  possible  to  use  data 
obtained  on  one  type  for  calculation  of  the  operation  of 
another  type.  A  coaplete  understanding  of  the  influence  of 
propellant  variables  on  drop  spray  has  not  yet  been  realized. 

An  experiaental  prograa  was  proposed  in  the  Seaiannual 
Report  (Ref,  2),  in  order  to  supply  needed  atoaization  data. 

The  objective  is  to  relate  propellant  properties  and  design 
paraaeters  to  droplet  spray  distribution.  Propellant  pro- 
perites  which  influence  droplet  distribution  are  density,  o,; 
viscosity,  yj;  and  surface  tension,  o.  Design  paraaeters  of 
importance  to  iapinging  atoaizers  are  injection  velocity,  Vj ; 
surrounding  gas  velocity,  u;  iapingeaent  angle,  0;  and  the 
injector  characteristics  (diaaeter,  D^;  length  ratio,  L/0^; 
and  outlet  configuration,  i,e,,  burrs  or  rounded). 
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S.  Soli4  Particl#  Effects  oa  Stability, 

T1i«  stabili zatioa  of  coabustioa  ia  solid  rockat  propellants 
coataiaiaf  powdarod  aotals  bas  stiaulatod  a  freat  deal  of 
research  ia  aa  effort  to  deteraiae  the  aechaaisas  of  the 
stabilizatioa  process.  Soae  of  the  aore  iaportant  observations 
coaceraiag  the  effects  of  aetal  powders  ia  solid  propellants 
aad  solid  exhaust  particles  are  suaaarized  in  Reference  2,  The 
principle  objective  of  this  review  of  solid  particle  effects 
(m  coabustioa  iastability  is  to  deteraine  possible  applica¬ 
tions  of  the  saae  principles  to  liquid  propellant  systeas. 

The  basic  equation  for  the  calculation  of  the  energy  dissipa¬ 
tion  for  conditions  applicable  to  liquid  systeas  is  given  by 
the  equation  of  Epstein  and  Carhart  (Ref.  16  ) . 


8 


6«crv 

— g —  (!♦*) 


_ i6z;i _ 

i6Z‘»472YZ*^«lY2(U22^2z2) 


The  actual  values  for  energy  lets  due  to  particle  attenua¬ 
tion  can  be  presi^nted  in  several  ways.  The  quantity  8  represents 

the  quantity  1  ;  which  corresponds  to  the  fractional  energy 

loss  per  foot  of  wave  travel*  Division  of  this  quantity  by  the 
wave  length  provides  the  energy  loss  per  cycle* 

In  an  aaalagous  aanner,  the  quantity,  Bxq  represents  E  \  ^dt 
or  the  fractional  energy  loss  per  second  and  division  by  the 
angular  frequency  represents  the  energy  loss  per  cycle* 


In  order  to  siaplify  the  presentation  of  results  the  quantity  3/ c 
is  presented*  This  quantity  represents  the  energy  loss  per 

E  article*  The  total  energy  loss  can  be  obtained  by  multiplication 
y  the  nuaber  of  particles  per  cubic  foot  of  volume.  For  a 
distribution  of  particles,  the  total  attentuation  can  be 
obtained  from  the  following: 


8  -  Z  8i 
1 


The  results  for  particles  corresponding  to  aluminum,  aluminum 
oxide,  hydrazine  and  mixtures  of  aluminua  and  hydrazine  are 
presented  in  figures  19  and  20  where  B/cw  is  plotted  as  a 
function  of  u  with  particle  radius  as  a  parameter.  When 
aultiplied  by.  q.  the  velocity  of  sound  in  the  gas.  the 
quantity  ^8/ cw  represents  the  energy  loss  due  to  attenuatior, 

I  dE.  n«r  particle  per  cycle*  It  is  evident  that  the 
E  'dt 

actual  attenuation  per  particle  increases  in  the  frequency 
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FIGURE  20.  VISCOUS  attenuation  PER  CYCLE  PETl  PARTI 


range  considered  as  the  particle  size  increases.  Particle 
density  also  has  an  effect,  the  heavier  particles  giving 
nore  attenuation.  This  is  illustrated  in  Figure  21  for 
10  nic^on  particles.  The  effect  of  particle  density  is 
appreciable  in  the  range  0  to  5000  and  becoaes  negligible  at 
frequencies  above  15,000. 

The  results  illustrated  in  Figures  19 , 20and  21  apply  for  a 
comparison  when  the  same  number  of  particles  of  each  size  or 
type  are  involved.  This  can  be  misleading,  however,  since  the 
factor  which  is  generally  held  constant  is  the  weight  of  solid 
materials.  For  a  given  weight  of  a  specified  material,  there 
are  more  small  particles  than  large  particles.  To  convert  the 
data  per  particle  to  data  per  unit  weight,  the  values  should  be 
divided  by  4/3  irr^.  In  Figure  22,data  for  aluminum  are  compared 
on  the  basis  of0/ccjr^  which  is  proportional  to  the  attenuation 
per  cycle  per  unit  weight.  Multiplication  of  the  factor  plotted 
in  Figure  22,  by  3  q  would  give  the  actual  attenuation  per 

T  »pp 

cycle  per  unit  weight  of  aluminum.  Above  a  frequency  of  2300,  5u 
particles  show  more  attenuation  than  lOw  particles  »imply  because 
there  are  8  times  as  many  Sy  particles. 

In  order  to  compare  the  attenuation  per  unit  weight  for  particles 
of  different  density,  the  density  terms  should  be  included.  A 
relative  comparison  is  given  in  Figure  23  where  (B/cwr^)  6  is 
plotted  against  frequency.  On  this  basis,  the  maximum  attenuation 
is  almost  independent  of  density  although  the  maximum  occurs  at 
different  frequencies. 

For  the  frequency  range  of  interest  in  most  liquid  systems, 
particle  sizes  of  the  order  of  2>20  microns  appear  to  be  of 
most  interest.  Although  particle  density  affects  the  actual 
value  of  the  energy  loss  due  to  particles,  the  differences 
are  not  great  so  that  hydrazine,  aluminum,  aluminum  oxide  and 
slurry  particles  have  similar  damping  effects.  * 


CYCLE  PER  PARTICLE,  4  Type  Particles,  10  Micron  Dia 


Aluminum 


FIGURE  22.  VISCOUS  ATTENUATION  PER  CYCLE  PER  UNIT  MASS. 


001 

693 

673 


111.  PARAHETRIC  STUDY  OF  COMBUSTION  INSTABILITY 

This  section  reports  the  results  of  a  paraeetric  study 
of  the  influence  of  liquid  rocket  engine  paraueters  on  coubustton 
instability.  A  basic  configuration  was  chosen  with  a  si»ilar 
geoKctry  as  the  Transtage.  Results  are  based  on  a  aono* 
actyAlyh/drasine  and  nitrogen  tetroxide  propellant  coabination* 
This  study  investigates  the  influences  of  propellant  aixture 
ratio»  injection  velocity,  diroplet  size  and  distribution,  and 
chaaber  pressure  on  the  ainiaua  pulse  strength  required  to 
trigger  instabi lit''. 

Results  of  this  section  are  obtained  by  a  aethod  using 
a  nonlinear  aodel  for  deteraining  the  zones  of  a  liquid  rocket 
engine  in  which  a  tangential  aode  of  high  frequency  instability 
is  aost  easily  initiated^  A  aethod  for  deteraining  these 
zones  was  developed  by  Beltran,  et  al  (Ref*  1  ).  A  rocket 
engine  is  analyzed  by  increaentally  dividing  the  coabustion 
chaaber  into  annular  nodes  in  the  r  and  z  directions*  Steady- 
state  properties  at  each  annular  node  or  position  in  the  chaaber 
are  coaputed  froa  the  propellant  vaporization  prograa*  These 
steady-state  properties  and  the  stability  liait  curves  froa 
the  instability  aodel  are  used  to  deteraine  the  stability  of 
the  node*  This  process  is  repeated  for  each  node  to  deter- 
Bine  a  stability  aap  of  the  entire  engine* 

1.  Steady-State  Analysis 

The  steady-state  perforaance  of  a  storable  propellant 
liquid  rocket  engine. of  Transtage  configuration  has  been 
investigated*  Injector  paraaeters  were  varied  in  a  coaputer 
paraaetric  study  which  utilized  the  Dynaaic  Science  Corporation 
steady-state  spray  coabustion  coaputer  prograa*  This  prograa 
(Ref.  1)  includes  a  dissociation  flaae  aodel  for  hydrazine  type 
fuels.  The  propellant  coabination  exaained  was  MNIl  and  N2O4. 

The  fuel  and  oxidizer  drop  sprays  are  characteristized  by  a 
logarithaiconoraal  distribution  of  drop  radii  about  a  specified 
aean  drop  size.  The  spray  is  divided  into  five  groups  of  equal 
mass.  The  group  of  saallest  radius,  however,  is  divided  into 
five  groups  again,  aaking  a  total  of  nine  groups  used  to  define 
the  spray.  This  was  done  to  obtain  wore  definition  in  the 
region  of  smallest  drop  size.  ■  The  mean  drop  size,  standard 
deviation  of  the  spray,  number  of  drop  groups,  and  mass  flow  rate 
are  the  parameters  used  to  define  the  spray  distribution. 

Table  1  shows  the  initial  drop  radii  and  number  of  drops  in 
each  of  the  groups  for  the  base  engine  spray  (spray  around  which 
the  injection  parameters  were  varied).  The  values  in  Table  I 
correspond  to  a  mass  mean  drop  radius  of  It  (  '03  in)  a  standard 
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deviation  of  2.3,  total  flow  rate  of  27.0  Ib/sec  and  a  2:1 
mixture  ratio. 


The  chamber  distance  [x]  along  which  the  governing  equations 
were  integrated  has  its  origin  at  the  atomization  point  of  the 
impinging  streams.  Figure  24  shows  a  water  flow  simulation 
of  the  stream  breakup.  The  impingement  diameter  of  the  streams 
is  .0625  inches  and  its  velocity  is  approximately  1200  in/sec. 
The  breakup  into  drops  is  seen  to  occur  within  O.SO  inches  of 
the  impingement  point.  Therefore,  the  assumption  of  an 
atomized  spray  vaporization  model  close  to  the  injector  face 
appears  justified. 


Tlie  parameters 
of  a  Transtage 
as  base  engine 


were  varied  around  values  selected  as  typical 
type  engine.  These  typical  values  will  be  denoted 
values.  (Shown  in  Table  11) 


Base  parameters  used  for  this  study  were: 


Propellant 

Injection  Velocity* 
Chamber  Pressure  (initial 
Initial  temperature* 

Total  Mass  Flow 
Mixture  Ratio 
Specific  Heat  Ratio  of 
Combustion  Gtses 
Mass  Mean  Drop  Radius* 
Standard  Deviation* 

Number  of  Groups* 


N2O4  and  MMH 
1000  in/sec 

value)  100  lb/in2 

530*  R 
27.0  Ib/sec 
2.0:1 

1.20 
.003  in 
2.3 
9 


The  chamber  configuration  was  specified  as  conical  with  a 
taper  from  107  in^  to  44  in^  in  a  17  in.  distance. 


A  series  of  computer  runs  were  made  with  specified  parameters 
varied  around  the  base  engine  values  shown  above,  bach 
parameter  was  varied  independently  of  the  others,  that  is, 
the  others  were  held  constant  at  base  engine  values.  This  was 
done  to  determine  the  significance  of  the  parameters  varied 
and  pinpoint  quantities  that  require  accurate  determination. 


The  first  two  runs  were 
The  runs  were  identical 
flame  was  suppressed  on 
vaporization  rate  [fract 


made  with  the  bb’^e  engine  quantities, 
with  the  excepti'n  that  the  dissociation 
the  second.  Tho  variation  of  the 
ion/inch]  of  tho  MMH  spray  with  distance 


Both  propellants. 


^  c  . 
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WATER  FLOW  SIMULATION  OF  STREAM  IMPINGEMENT 


FIGURE  24.  Water  Streaa  Inpinging  at  90*. 

Stream  Diameter  .062S  inch 
Injection  Velocity  Approximately  1200 


in/scc. 
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TABLE  II 

STEADY-STATE  «ad  INSTABILITY  PABANETERS 


SI 

Varied 

Paraaeters 

CoBtroliinf 

Propellaat 

NiaiauB 

AV  laches 

Nils 

L 

AP/P 

1 

Base  Eufine 

OX 

0.24 

0.94 

320 

0.40 

0.032 

D 

Nithout  TF 

Fuel 

0.7S 

0.75 

379 

0.19 

0.050 

H 

2Sv 

OX 

0.03 

0.30 

105 

4.62 

0.080 

■ 

22Sw 

OX 

1.43 

3.04 

1021 

0.05 

0.099 

H 

lo 

OX 

0.76 

2.95 

1234 

0.10 

0.C37 

B 

lo 

OX 

0.03 

0.24 

93 

3.24 

0.087 

B 

300  psi 

OX 

0.75 

1.27 

1361 

0.12 

0.032 

8 

SOO  psi 

Fuel 

1.53 

1.68 

3270 

0.05 

0.052 

9 

SOO  in/sec 

OX 

0.07 

0.88 

305 

0.75 

0.032 

10 

2000  in/sec 

OX 

0.89 

1.01 

352 

0.23 

0.040 

11 

0/F  1.6/1 

OX 

0.28 

0.95 

431 

0.31 

0.030 

12 

0/F  2.5/1 

OX 

0.25 

0.93 

316 

0.48 

- 

0.031 

<•««  tkt  chMk«r  froB  tli«  pviat  of  otoBisotioii  is  skowB  in 
Fifvro  2S.  loelofiog  tko  oxotkonic  4isseciotio«  fUo«  produces 
•pproxiMtolp  •  SOO-400%  iocrooso  la  tko  iaitial  vaporisatioa 
rate.  After  2.S  iaekos  tko  rates  are  seea  te  cress  as  tke  kifher 
iaitial  rate  kas  caused  tke  less  ef  tke  sualler  drops  frou  tke 
distrikutiea  leaviag  ealy  tke  sleuer  vaporisiag  larger  sise  drops. 

At  tke  cross-over  peiat  (x«2.S  iackes),  tko  iaclusioa  of  tke  two 
flaae  iwdel  kas  raised  tke  total  fractioa  of  tke  fuel  spray  vaporized 
frea  20%  to  42%. 

Tko  varii^tioa  of  tko  oxidisor  vaporization  rate  witk  chaaber 
distance  for  tkoso  two  runs  is  skown  ia  Figure  26.  Tke  change 
is  ratkor  slight,  however,  there  is  sene  variation  since: 

a.  The  higher  fuel  vaporisation  rate  causes  the  chaaber 
gas  velocity  to  increase  uore  rapidly  affecting  the  heat  and 
Bass  transfer  coefficients  to  the  drop. 

b.  Tke  local  0/F  ratio  is  altered  by  the  high  fuel 
vaporisation  rate.  This  affects  the  adiabatic  flaae  tenpera- 
tore  (driving  potential  for  heat  transfer]  as  well  as  the 
partial  pressures  of  the  species  in  tke  conbustion  gases 
(driving  potential  for  aass  trvVnsfer]. 

c.  Tke  change  ia  local  0/F  ratio  between  these  runs 

also  causes  changes  in  the  transport  properties  of  the  diffusion 
aantle  by  altering  the  equilibriua  conposition  of  the  conbustion 
gases.  Note:  Thus  the  fuel  vaporisation  rate  is  increased 
significantly  by  the  two  flaae  aodel  and  this  effect,  in  turn, 
causes  an  increase  in  the  oxidiser  vaporization  rate;  the  two 
being  interrelated.  Of  particular  interest  is  the  change  in 
the  local  0/F  ratio  with  the  inclusion  of  the  two  f lanes.  The 
calculat  ^  values  of  0/F  are  shown  in  Figure  27  for  these  two 
runs.  The  two  flane  aodel  has  caused  the  conbustion  process 
to  switch  froB  fuel  controlling  (0/F>2.5)  to  oxidizer  controlling 
[0/F<2.S].  For  reference,  the  nachine  plots  of  fraction 
vaporised,  gas  Mach  nuaber,  fuel  and  oxidizer  drops  radii, 
and  fuel  and  oxidizer  drop  velocities  for  the  runs  with  and 
without  the  two  flaae  burning  regiae  are  shown  in  Figure  28-35 
and  34-39  respectively.  It  was  noted  that  the  fuel  drops  were 
at  all  tiaes  within  the  two  flaae  regiae  for  the  base  engine 
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FIGURE  27.  Effect  of  Fuel  Dissociation  Flaae  on  Local  0/F  Ratio 
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PIptiRE  38.  Pud  Drop  Velocities  (n©  doconposition  flatne) 
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FIGURE  :9.  Qxldiier  Drop  Velocities  (using  decomposition  flame  model) 
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Chanb^r  pressure  was  varied  fron  the  base  value  of  100  psia 
to  300  and  SQO  psia.  The  effects  on  the  ^uel  and  oxidizer 
vaporization  rates  are  shown  in  Figures  40  and  41.  Increasing 
the  chatabor  pressure  fron  100  to  300  psia  increases  the  vapori¬ 
zation  rate  of  the  fuel  and  the  oxidizer  with  the  greatest 
effect  being  felt  by  the  oxidizer.  This  can  be  seen  in  the 
plot  of  0/F  ratio  shown  in  Figure  42.  Through  the  sensitive 
zone  (zero  relative  velocity)  the  0/F  ratio  at  100  psia  is  about 
1.0  while  at  300  psia  it  is  l.S.  Raising  the  ehanber  projsuro 
to  500  psia  causes  the  oxidizer  vaporization  rate  to  further 
increase.  At  500  psia  the  boundary  layer  thickness  over  tho 
throe  snallest  drop  groups  is  sufficiently  snail  so  that  the- 
oxidation  and  dissociation  flacos  occur  at  tho  sano  position* 
This  noans  the  fastest  vaporizing  drops  aro  out  of  tho  two 
flano  rogine  and  will  vaporize  correspondingly  slower*  In 
Figure  40  the  fuel  vaporization  rate  at  500  psia  is  soon  to 
bo  lower  than  at  100  and  300  psia  until  tho  snallest  drop 
group  [tho  group  controlling  the  initial  vaporization]  is 
consuned.  At  this  distance  (xa0.*35]  the  drops  aro  in  tho  two 
flano  regino  since  the  boundary  layer  thickness  approaches 
infinity  as  tho  relative  velocity  approaches  zero.  Thcreforot 
the  fuel  vaporization  rate  undergoes  alnost  a  stop  chango» 
decreasing  rapidly  at  xa0.40,  as  the  controlling  drop  group 
is  consuned*  Tho  0/F  ratio  undergoes  largo  excursions  in 
this  region  as  can  be  seen  in  Figure  42*  Of  particular  interest 
is  tho  fact  that  tho  fuel  has  bocone  tho  controlling  propellant 
at  tho  sensitive  zone.  Thoroforo»  increasing  oporating  pressure 
of  tho  engine  can  switch  control  of  tho  conbustion  process  to 
tho  fuel. 

Varying  tho  nass  noan  drop  size  has  a  largo  offoct  on  the 
vaporization  rate  of  both  the  fuel  and  tho  oxidizor*  Figures 
43  and  44  show  the  fuel  and  oxidizer  vaporization  rates  for 
mass  noan  drop  radii  of  25,  75,  and  225  sprays*  Tho  standard 
deviation  of  tho  distribution  was  naintainod  at  2.3.  Tho  drop 
radii  and  nunbor  in  tho  groups  describing  tho  sprays  aro  shown 
in  Table  HI  for  tho  25u  and  225{j  sprays*  Tho  fraction  vapori¬ 
zed  of  tho  fuel  and  oxidizor  sprays  aro  shown  in  Figures  45  and 
46.  The  vaporization  rate  of  tho  25u  spray  is  soon  to  bo  nuch 
faster  than  both  tho  7Sv  and  225p  sprays  because  of  tho  largo 
nunbors  of  snail  drops*  For  roforonco  tho  variation  in  tho 
drop  radii  conprising  tho  fuel  spray  distributions  is  shown 
in  Figures  47  and  48  for  tho  2S)i  and  225u  fuel  sprays*  Tho 
equivalent  plots  for  tho  oxidizor  sprays  aro  shown  in  Figures 
49  and  50.  In  all  throe  cases  tho  0/F  ratio  plot  [Figure  51] 
shows  tho  oxidizor  to  bo  controlling  tho  conbustion  process* 


nation  Rate  with  Chamber  Press 


ICURE  41^  Effect  of 


Chaober  Distance  [x] 
FIGURE  42.  Variation  of  0/F  Ratio  with 


TABLE  III 


OXIDIZER  SPRAY  DISTRIBUTION 
o«2.3,  00*18.  Ib/sec. 


-  2Su 

(.001  in) 

Tq  b  22Sti 

(.009  in) 

Group 

Radius 

(oil)* 

Nunber 

Radius 

(ail)  ' 

Nuober 

1 

0.1807 

4.430x1011 

1.6259 

6.077x10® 

2 

0.2740 

1.269x1011 

2.4664 

1.741x10® 

3 

0.3443 

6.399x101® 

3.098? 

8.778x10^ 

4 

0.4073 

3.865x101® 

3.6659 

5,302x10^ 

S 

0.4675 

2.S57xl0l® 

4.2072 

3.507x10’ 

1 

6 

0.6476 

4.808xl0l® 

5.8285 

6.596x10’ 

7 

1.0000 

1.306x101® 

9.0000 

1.791x10’ 

8 

1.5441 

3.547x10® 

13.8973 

4.866x10® 

9 

2.9044 

5.330x10® 

26.1397 

7,312x10® 

FUEL  SPRAY  DISTRIBUTION 
0*2.3,  nf*9.0  lb/5cc. 


*  2Sw 

(.001  in) 

Tjj  »  225y 

(.009  in,) 

Group 

Radius 

(nil) 

Number 

Radius 

(nil) 

Number 

1 

0.1607 

2.805x1011 

1.6259 

3.847x108 

2 

0.2746 

8.034x101® 

2.4664 

1.102x108 

3 

0.3443 

4.051x101® 

3.0987 

5.557x107 

4 

0.4073 

2.447x101® 

3.6659 

3.556x10? 

5 

0.4673 

1.6i9xl0l® 

4.2072 

2.220x10? 

6 

0.6476 

3.044x101® 

5.8285 

4.176x10? 

7 

1.0000 

8.268x10® 

1.134x10? 

8 

1.5441 

2.246x10® 

13.8973 

3.080x10® 

9 

BOB 

26.1397 

4.629x10® 

*1 
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Ti^(3  effects  ef  initial  drop  volocity  on  tho  fuel  and  oxidizer 
sprey  vaporization  rate  are  shoun  in  Figures  S2  and  S3.  Tho 
injoetion  volocitios  of  tho  sprays  wore  varied  sinultanoously 
around  tho  baso  value  of  lOt)  in/soc»  Tho  selected  values 
were  500  and  2000  in/seo«  Even  at  a  2000  in/sec.  injection 
velocity  all  drops  of  tho  fuel  distribution  were  in  the  two 
fleoo  burning  rogin*-  increasing  tho  initial  velocity  causes 
tho  vaporization  rato  to  decrease  as  tho  drop  dwell  tine  in 
each  sognont  of  the  chacber  is  shorter.  The  nachino  plots 
of  tho  fraction  of  the  fuol«  oxidizer,  and  total  spray 
vaporized  is  shown  in  Figures  S4,  5S,  and  S6  for  injectioa 
volocitios  of  SOO,  1000,  and  200Q  inches/scc.  As  an  exanplo 
the  oxidizor  fraction  vaporized  at  1.0  inch  is  .30,  .20,  and 
.12  for  SOO,  1000,  and  2000  inch/soc.  respectively.  The 
equivalent  values  for  tho  fuel  spray  are  .34,  .27,  and  .18. 

Tho  longer  dwell  tine  in  tho  chanbor  distance  of  1.0  inches 
with. tho  SOO  inch/soc.  injection  velocity  is  reflected  by 
tho  larger  vaporized  fraction.  Another  effect  of  varia¬ 

tion  of  tho  injection  velocity  is  tho  distribution  of  the  zero 
rplativo  velocity  points  as  the  injection  velocity  is  increased. 
Tho  nachino  plots  of  tho  drop  velocities  with  gas  velocity 
superieposod  are  shown  in  Figures  S7,  58,  and  59  for  tho  oxi¬ 
dizer  sprays  injected  at  SQO,  1000,  and  2000  inch/sec*  The 
corresponding  plots  for  tho  fuel  sprays  are  shown  in  Figures 
60,  61,  and  62.  Tho  gas  velocity  increases  sufficiently  fast 
so  that  all  drops  go  through  tho  zero  relative  velocity  point 
at  approxioatoly  O.lO  inches  when  injected  at  500  inchos/soc. 

At  this  velocity  drag  i»  relatively  low  and  the  drops  go  through 
the  zero  point  [Figures  57  a.,  i  60)  before  slowing  down.  With 
un  injection  velocity  of  1000  inches/socond  tho  drops  go  through 
the  zero  point  [Figures  58  and  61]  at  0,20  for  the  smallest 
group  and  0,25  inches  for  tho  largest  drops.  The  distribution 
still  goes  through  the  zero  point  at  relatively  tho  same  point. 

At  a  2000  inchos/soc.  injection  velocity,  however.  Figures  59 
end  62  show  tho  snallost  drops  to  go  through  the  sensitive  region 
at  0.60  inches  while  tho  largest  reaches  its  zero  relative' 
velocity  point  at  1.0  inches.  This  effect  of  a  spreading  of 
tho  sensitive  region  complicates  tho  calculation  of  combustion 
stability  by  a  single  drop  model  and  necessitates  the  deter¬ 
mination  of  an  effective  drop  that  characterizes  the  spray. 

In  tho  next  sot  of  runs  tho  standard  deviation  of  tho  drop 
radii  distribution  around  the  75u  moan  was  varied  from  the 
base  value  of  2.3  to  values  of  4.0  and  1.0.  A  value  of  1*0 
pjroscribos  a  spray  of  uniform  drops.  A  summary  of  the  drop 
radii  and  number  of  tho  groups  describing  this  spray  is  shown 
in  Toblo  IV.  Onco  again  tho  first  five  groups  of  the  distri¬ 
bution  represents  20%  of  tho  mass.  Tho  offoct  of  standard 
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OXIDIZER  SPRAY  DISTRIBUTION 
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deviation  on  the  vaporization  rate  of  the  total  fuel  and 
oxidizer  spray  is  shown  in  Figures  63  and  64,  hith  a  standard 
deviation  of  4.0  the  radii  distribution  is  spread  out  making 
"lore  small  drops.  Therefore*  the  initial  vaporization  rate 
is  much  higher  for  this  case.  ,As  the  smaller  drops  are  con¬ 
sumed*  however*  the  rates  approach  each  other  since  the  larger 
drops  become  controlling.  The  fuel  vaporization  rate  with  the 
deviation  of  1.0  is  relatively  constant  as  can  be  seen  in 
Figure  63.  The  calculated  variation  of  the  local  0/F  with 
distance  down  the  chamber  is  shown  in  Figure  65.  It  is  seen 
that  the  oxidizer  is  controlling  for  approximately  4,0  inches. 
Large  differences  in  the  spray  fraction  vaporized  are  noted 
between  these  runs  out  to  a  chamber  distance  of  about  2  inches. 
After  this  length  the  standard  deviation  has  little  effect  on 
the  fraction  vaporized. 
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FIGURE  64,  Variation  of  Fuel  Vaporization  Rate  with  Standard  Deviation  of  Spray. 
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FIGURE  6S.  Variation  of  0/F  Ratio  with  Standard  Deviation  of  Spray 


2.  Stability  Analysis 

Coabustion  instability  results  are  obtained  by  a  aethod 
using  a  nonlinear  aodel  for  deteraining  the  zones  of  a  liquid 
rocket  engine  in  which  a  tangential  node  of  high  frequency  in¬ 
stability  is  aost  easily  initiat3d«  The  sbthod  for  deteraining 
these  zones  was  developed  by  Oynaeic  Science  Corporation  and 
details  are  presented  in  Reference  1.  This  aethod  uses  a  non¬ 
linear  instability  aodel  with  a  propellant  vaporization  prograa. 

A  rocket  engine  is  analyzed  by  increnentally  dividing  the  con- 
bustion  chanber  into  annular  nodes  in  the  r  and  z  directions  as 
illustrated  in  Figure  &•  Steady-state  properties  at  each  annular 
node  or  position  in  the  chanber  are  conputed  froa  the  propellant 
vaporization  program.  These' steady-state  properties  and  the 
stability  limit  curves  fron  the  instability  model  are  used  to 
determine  the  stability  of  the  node*  This  process  is  repeated 
for  each  node  to  determine  a  stability  map  of  the  entire  engine. 
Details  of  the  application  of  this  method  are  presented  in 
Reference  1.  This  section  will  concentrate  on  the  results  of  a 
parametric  study  of  the  influence  of  liquid  rocket  engine  para¬ 
meters  on  combustion  instability.  This  study  investigates  the 
influences  of  propellant  mixture  ratio,  injection  velocity, 
droplet  size  and  distribution,  and  chamber  pressure  on  the 
minimum  pulse  strength  required  to  trigger  instability. 

Steady-state  and  instability  parameters  at  the  most  sensi¬ 
tive  region  in  the  engine  are  presented  in  Table  jt.  The  para¬ 
meters  varied  were  chosen  as  an  average  of  what  night  bo  expoctod 
on  a  "Transtage"  type  engine,  rather  than  to  represent  exact 
values  actually  obtained.  These  parametric  variations  around 
the  base  configuration  were: 

0/F  1.6  and  2.5 

300  and  500  psi 

Mass  median  droplet 
radius  25u  and  225v 

Droplet  distri¬ 
bution  1  and  4o 

Injection  velo¬ 
city  iOO  and  2000  in/sec. 

Results  of  the  parametric  study  show  that,  for  the  conditions 
considered,  the  oxidizer  controlled  the  combustion  process 
except  where  the  "two  flame"  nodol  was  olininatod  or  the 
chanber  pressure  was  above  300  psi*  The  "two  flames"  wore 
suppressed  in  the  500  psi  case  due  to  the  thin  convective  filn 
thickness  at  high  chamber  densities.  The  study  shows  that  the 
propellant  vapor  phaso  reactions  reverse  the  vaporization  rate 
trend  indicated  by  the  heat  of  vaporization,  i.o.  the  ozidizor 
vaporizes  slower  than  the  fuel,  thus  enabling  the  oxidizer  droplets 
to  control  the  combustion  process. 
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^QS'alts  show  thst  undor  the  conditions  studied  droplet 
si so  end  distribution  hsvo  the  largest  effect  on  stability. 

Figu/ros  dd  end  5?  present  the  variation  of  niniaun  pulse 
s^trongth  versus  oass  sean  droplet  radius  and  distribution. 

For  a  cheober  configuration  of  the  Transtage  there  is  an 
optimum  Bass  codian  droplet  size  for  instability  around  75ii. 
Incroesing  or  doeroasing  the  droplet  size  increases  the  sta« 
bility.  For  very  largo  droplet  radii  the  burning  rate  paraceter 
is  very  soall  end  the  droplet  Reynolds  nunber  is  large;  however, 
the  not  effect  is  to  fall  to  tho  ^aft  side  of  the  &F/P  vs.  L 
euryo  rosuiting  in  incroasod  stability.  For  very  snail  droplet 
radii  tho  burning  rate  paranetor  is  large;  however,  tho  Reynolds 
nunbor  is  low,  thus  any  burning  rate  parasoter  will  have  a  high 
dP/F.  Thoro  is  also  an  optiaun  distribution  for  instability  at 
approzioatoly  2a,  increasing  the  distribution  to  4o  and  decreasing 
tho  distribution  to  lo  increas«»4  stability.  Increases  in  stability 
occur  for  tho  sans  reasons  as  for  the  effects  of  droplet  radius. 

The  greatest  stabilizing  effect  was  obtained  at  the  large  dis¬ 
tribution  end.  At  4o  thoro  aro  a  great  nunber  of  very  snail 
and  largo  droplets,  with  tho  snail  ond  controlling  tho  stability 
charaetor  of  tho  spray.  Snail  droplets  produce  nore  stable 
coobustion  since,  while  tho  stoady-stato  burning  rate  is  high, 
thoir  burning  response  w'  is  low,  thus  decreasing  the  wavo- 
conbustion  rate  coupling. 

Increasing  tho  droplet  size  noved  the  nininun  relative 
velocity  region  further  fron  tho  injector  face,  while  decreasing 
tho  droplet  distribution  hsd  tho  sano  effect.  This  is  just  the 
result  of  incroasod  droplet  surface  area  duo  to  a  greater  nunber 
of  snallor  droplets  for  a  given  flow  rate.  It  is  interesting  to 
note  that  a  cass  nod^  !!in  droplet  sizo  of  7Sw  and  a  distribution  of 
2.3o  correspond  to  tho  worst  conditions  for  stability  and  also 
droplet  sizo  and  distribution  pToducod  by  typical  Transtage  type 
injectors.  Figure  68  has  boon  dovolopod  fron  tho  results  of 
Roforonco  3  with  corrections  for  MMH  and  NTO  physical  properties. 
This  figure  shows  that  parallel  and  inpinging  jets  have  distribu¬ 
tions  of  2.3o,  while  in-lino  triplets  have  a  distribution  of  3.6o. 
For  tho  study  conducted,  i.o.,  base  engine  conditions  with  a  75p 
nass-nodian  droplet  sizo,  a  triplet  injector  would  bo  norc  stable 
than  parallel  or  doublet  jots.  Fron  Figure  68,  droplet  sprays 
botwoon  SOw  to  lOOu  aro  produced  by  jot  dianotors  for  parallel 
jots  (-0.00S  to  0.025  inches),  inpinging  jots  (-0.035  to  0.095 
inches),  and  triplets  (-0.125  to  0.650  inches).  Injector  atoni- 
zatioh  photographs  taken  by  Dynanic  Science  Corporation  indicate 
that  A  positioned  triplets  produce  a  coarser  spray  that  in-lino 
triplets  and  X  quadlots  produce  a  coarsor  spray  than  A  triplets. 
This  is  duo  to  tho  fact  tho  doublets  and  in- lino  triplots  produce 
fans.  Those  fens  fern  droplets  by  tho  fornation  of  shoot  insta¬ 
bilities  with  finer  droplets  being  produced  by  thinner  shoots. 
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Droplet  Distribution  V 

FIGURE  67.  Efxoct  of  Droplot  Distribution  on  Mininun  Threshold  Disturbance 
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Quadlots  end  a  triplets  tend  to  fors  a  stroan  which  broahs 
up  siailar  to  parallel  jets  thus  foroing  droplet  aizos  between 
the  in-lino  triplet  and  parallel  jets. 

Chacber  pressure  was  the  third  nost  important  paraooter 
affecting  instability.  Increasing  chaober  pressure  aoved  the 
siniouQ  relative  velocity  position  away  frpn  the  injector 
face.  This  occurs  due  to  tho  fact  that  the  gas  velocity 
increases  slower  than  the  increased  nass  addition,  since  the 
gas  density  is  higher.  At  SOO  psi  tho  burning  rate  was  slower 
since  tho  MMH  droplets  burned  with  a  single  flane.  At  this 
pressure  the  gas  density  is  high  thus  the  gas  Reynolds  nunber  is 
high,  resulting  in  a  thin  diffusion  fila  thickness.  Since  tho 
diffusion  tine  is  short  due  to  the  thin  fila  tho  MMH  does  not 
havo  ties  to  decoapose.  Thus  tho  heat  transfer  is  lower  to  the 
droplet  and  the  vaporization  process  becoaos  slower.  Tho  effect 
of  chaebor  pressure  on  stability  is  shown  in  Figure  69.  For 
constant  flow  rates,  i.e.,  the  contraction  ratio  is  increased  to 
increase  chaaber  pressure,  stability  roaains  constant  with  in¬ 
creased  chaabor  pressure  between  100  and  300  psia.  Hero  tho 
stabilizing  effect  of  decreasing  the  burning  rate  paranetor  by 
increasing  the  contraction  ratio  equals  the  destabilizing  effect 
of  increasing  tho  Reynolds  nuaber  of  the  droplet  by  increasing 
chaabor  gas  density.  For  the  conditions  considered,  increasing 
the  chaabor  pressure  (at  constant  flow  rate)  to  500  psia  pro¬ 
duced  a  stabilizing  effect.  Since  at  this  pressure  the  droplet 
spray  considered  burns  in  tho  slower  single  flane  regi'no.  Thus, 
the  fraction  vaporization  rate  decreases,  decreasing  the  burning 
'rate  paranetor  offsetting  tho  destabilizing  effect  of  the  gas 
density.  For  larger  droplet  sprays,  whore  the  droplets  would 
continue  to  burn  in  the  two  flane  regine  at  500  psia,,  or  pro¬ 
pellants  which  do  not  decoapose,  stability  would  ronain  constant 
or  decrease.  Tho  500  psia  case  was  tho  only  condition  where 
control  switched  froa  the  oxidizer  spray  to  the  fuel  spray. 

To  dotoraine  tho  differences  between  MMH  burning  with  a  "two 
flane"  nodol  and  a  single  flaae  considered  in  Roference  3,  a 
case  was  run  at  base  conditions  with  the  deconposition  flane 
suppressed.  For  this  run  tho  burning  rate  paranetor  was  reduced 
froQ  0.40  to  Q.IO  corresponding  to  an  increase  of  ainiaun  threshold 
disturbance  fron  0.032  psi  to  .050  psi.  Thus  additives  which 
would  suppress  the  MMH  docoaposition  would  increase  stability; 
howovor,  this  nay  tend  to  decrease  perfornanco. 

For  conditions  whore  tho  chaabor  pressure  is  increased 
by  increasing  tho  flow<  rate,  while  keeping  tho  contraction 
ratio  constant,  tho  engine  is  dostabilizod.  Under  theso  con¬ 
ditions  tho  increased  Reynolds  nuaber  of  tho  droplet  overrides 
the  decrease  of  the  burning  rate  paranetor  due  to  tho  decreased 
fraction  vaporized  per  inch. 
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In^octioQ  velocity  effects  are  shovn  in  Figure  70. 

The  point  of  ninieun  relative  velocity  was  coved  away  free 
the  injector  face  by  increasing  the  injection  velocity. 
Increasing  the  injection  velocity  decreases  the  stay  tico  of 
the  droplets  requiring  core  distance  to  vaporize  the  sane 
acount.  As  shown  in  Figure  70  increasing  the  injection 
velocity  slightly  increases  the  stability.  The  increased 
injection  velocity*  decreases  the  fraction  vaporised  per  inch 
and  thus  the  burning  rate  paracotor, resulting  in  increased 
stability. 

The  effects  of  eixturo  ratio  are  presented  in  Figure  71. 
There  is  essentially  no  effect  of  nixturo  ratio  on  disturbance 
aoplitudo  which  rooains  approxisately  at  0.030  psi.  It  should 
bo  noted  that  this  variance  of  nixturo  ratio  only  represents 
the  proportions  of  propellant  with  the  total  flow  rate  kept 
constant o  Exporinontally ^strong  effects  of  nixturo  ratio  are 
noted  on  stability  characteristics i  however*  when  nixturo  ratio 
is  veviod  in  a  given  configuration*  changes  result  in  droplet 
spray*  injection  velocity  and  chanbor  pressure  as  well  as  local 
nixturo  ratio.  Thus*tho  effects  of  changing  nixturo  ratio  in 
an  engine  result  free  the  influence  on  other  paranetors* 

Expcrinontal  Correlation 


Rocket  engine  tests  utilising  the  MNH/nitrogen  totroxido 
propellant  conbination  were  conducted  at  the  Air  Force  Rocket 
Prepulsion  Laboratory  and  results  wore  reported  in  Roforonco  17.. 
Twenty  throe  tests  were  conducted  using  two  Transtago  engine 
configurations.  The  two  engines  consisted  of  the  Titan  Ill 
Transtago  coebustion  syston*  i.o.*  thrust  chanbor*  injector 
and  bipropollant  valve.  The  stability  evaluation  tost  progran 
used  a  workborso  stool  thrust  chanbor  having  the  sane  conical 
internal  contour  as  the  flight  nodol  ablative  chanbor.  Models 
Trax  21- ID  and  Trax  21-llB  injectors  wore  tested.  Major  effort 
was  oxpondod  on  the  Trax  21-lD  injoctor  which  was  tested  14  tines 
with  at  chanbor  pressures  between  91-102  psia  and  nixturo 
ratios  between  1.75-2.56.  The  Trax  21-llB  was  tested  4  tines 
with  MtCH  at  ehanber  pressures  between  97-100  psia  and  nixturo 
ratios  between  1.78-2.33. 

As  described  in  Reference  17*,  the  -ID  injoctor  has  an 
aluninun  dished  face*  is  11.645  inches  in  dianotor*  and  possesses 
four  radial  baffle  vanes  extending  fron  a  snail  dianotor  center 
fuel  hub.  All  the  fuel  is  used  to  rogenorativoly  cool  the  baffles 
prior  to  injection  into  the  conbustion  ehanber.  The  -ID  injoctor 
has  a  quadlot  pattern  consisting  of  two  oxidizer  holes  inpinging 
bitb  two  fuel  holes.  The  total  included  inpingonont  angles  for 
both  sots  of  holes  is  70  degrees.  There  are  a  total  of  8  rows 
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of  injector  eloeonts  with  three  fuel  holes  in  each  baffle  tip 
and  one  row  of  fuel  showerhead  holes  located  around  the  peri¬ 
phery  of  the  injector  provide  for  filn  cooling.  There  aTO 
672  oxidizer  orifices  with  an  orifice  diacotor  of  0.0360  inches 
and  672  fuel  orifices  with  an  orifice  diaooter  of  0.0292  inches. 

The  Trax  21-llB  injector  pattern  also  consists  of  8  row^ 
of  injector  elecents  of  which  the  first*  second,  third,  fifth, 
and  seventh  rows,  starting  fron  the  injector  center,  are  the 
quadlet  pattern  as  the  -10  injector.  The  fourth,  sixth,  and 
eighth  rows  were  changed  to  provide  for  A  triplets  with  two 
fuel  holes  ispinging  with  one  oxidizer  hole.  The  fuel  hole 
'diasotor  has  been  enlarged  to  0.0512  inches  (0.0360  inches  for 
the  quadlet  pattern).  In  addition,  tho  outer  rows  of  fuel  fila 
coolant  holes  have  been  eliminated. 

The  -10  and  -118  injectors  had  approximately  the  same 
performance,  with  the  -118  being  higher  at  a  mixture  ratio  of 
1.8.  All  14  tests  conducted  with  tho  Transtago  -D  injector 
resulted  in  combustion  instability;  three  of  these  instabilities 
were  initiated  by  spontaneous  pops  with  a  AP  higher  than  40  psi, 
and  the  remainder  of  the  tests  were  initiated  with  10,  15,  or 
20  -grain  pulse  charges.  The  4  tests  conducted  with  the  -118 
injector  also  wore  triggered  unstable. Throe  of  tho  4  tests 
required  a  40-grain  pulse  to  initiate  instability  while  tho 
fourth  required  only  a  20-grain  pulse.  Pulse  charge  size 
required  to  trigger  instability  for  the  Transtago  injectors 
is  plotted  in  Figure  72.  All  resultant  instabilities  with 
tho  Transtago  injectors  wore  of  tho  second  tangential  mode  with 
frequencies  of  oscillation  of  approximately  3700-3900  cps. 
Peak-to-poak  pressure  amplitudes  wore  higher  for  the  -ID 
injector  (100-130psi)  than  the  -118  injector  (40-60  psi). 

In  general,  the  11-8  injector  showed  much  greater  stability 
than  the  -10  injector. 

To  determine  the  droplet  sizes  produced  by  tho  -10  dnd 
-118  injectors,  water  flow  tests  were  conducted  at  Dynamic 
Science  Corporation.  By  taking  spray  photographs  it  was 
determined  that  X  quadlets  produced  larger  droplets  than 
doublet  injectors  and  A  triplets  produced  larger  droplets 
than  in-line  triplets  under  tho  samo  conditions.  Larger 
droplets  are  formed  si.ico  quadlets  tend  to  produce  a  stream 
rather  than  a  thin  sheet,  resulting  in  wave  instabilities  of 
longer  period  on  the  liquid  surface.  In-lino  triplets 
form  fans,  while  A  triplets  tend  to  form  thicker  shoots  or 
stream's.  Using  tho  photographs  taken  by  Dynamic  Science  and 
results  of  Figure  68 ^mass-nodian  droplet  sizes  of  the  quadlet 
and  triplet  are  estimated  to  bo  65v  and  20u  rospoctivoly . 


The  raass-nedian  fuel  droplets  produced  by  the  cooling  fils 
parallel  jets  was  estisated  to  be  llSy.  Since  the  outer  edge 
of  the  engine  is  nest  sensitive  to  be  pulsed  unstable,  the 
effective  droplet  aass-nedian  size  can  be  estisated  to  be 
approxinately  70u  for  the  *10  injector  and  30u  for  the  -IIB 
injector  by  weighting  the  outer  elenonts  heavier  than  the 
center.  Pros  Figure  66  it  can  be  seen  that  for  the  propellant 
and  conditions  of  the  Transtage  engine  the  -IIB  injector  would 
be  noro  stable.  In  addition,  the  >11B  injector  also  has  a 
greater  droplet  distribution  since  triplets  are  being  utilized. 
Greater  droplet  distribution  improves  stability  for  the  case 
considered.  As  shown  in  Figure  72  there  is  a  scatter  in  the 
minimus  pulse  strength  required  to  produce  instability.  Since 
these  runs  were  made  at  diiferent  mixture  ratios,  variation 
occurred  in  injection  velocity,  chamber  pressure,  and  local 
mixture  ratio.  Secondary  effects  also  occur  on  the  atomization 
process.  Thus  for  these  two  injectors  considered,  injector 
patterns  influenced  the  combustion  characteristics  considerable 
more  than  operating  conditions. 

Conclusions 


Before  a  detailed  experimental  correlation  can  be  obtained 
greater  information  about  injector  spray  characteristics  will  bo 
required.  While  correlation  was  obtained  with  the  two  injector 
tests  considerably  more  work  will  be  required  to  conclusively 
evaluate  the  value  of  the  outlined  model.  While  trends  seem 
to  be  predicted  there  is  a  largo  difference  between  predicted  and 
measured  AP/P  required  to  trigger  instability.  It  was  experi¬ 
mentally  determined  that  the  pulse  amplitude  required  to  trigger 
instability  was  0.4  while  calculated  amplitudes  are  approxinately 
0.03.  Greater  refinement  of  the  model  will  be  required  before 
quantitative  results  can  bo  obtained;  however  a  method  which 
gives  qualitative  results  is  of  value  from  a  design  standpoint. 

Significant  results  of  this  study  show  that:  1)  the  two 
flame  model  showed  increased  burning  rate  and  decreased  stability. 

The  oxidizer  controlled  the  combustion  process,  except  at  500  psi 
where  the  two  flames  were  suppressed,  thus  the  propellant  vapor 
phase  reactions  reverse  the  vaporization  rate  trend  indicated  by 
the  heat  of  vaporization.  2)  There  was  a  droplet  size  and  distribu¬ 
tion  which  produced  minimum  stability.  Variation  in  either 
direction  produced  increased  stability.  3)  Chamber  pressure 
variations  resulting  from  changes  in  contraction  ratio  at  constant 
flow  rate,  increased  the  stability.  Keeping  the  contraction  ratio 
constant  and  increasing  chamber  pressure  by  increasing  flow  rate 
decreased  stability.  4)  Stability  increased  with  increasing  injection 
velocity.  Spray  characteristics  were  held  constant  while  injection 
velocity  was  varied.  5)  Variance  of  mixture  ratio,  resulting  from 
changing  the  proportions  of  propellant  with  constant  flow  rate, 
did  not  affect  the  stability  characteristics.  6)  Comparison  with 
experiaents  showed  that  the  model  underestimated  minimum  disturbance 
levels  but  correctly  indicated  stability  trends  with  injector  designs. 
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IV,  LIQUID  PHASE  MIXING  AND  REACTION 


Tho  acount  of  reaction  between  inpinging  unlike  streams 
oust  bo  considered  when  calculating  the  combustion  rate  of  hyper- 
golic  propellant  systems.  Liquid  phase  reaction  and  heat  release 
during  impingement  would  affect  the  droplet  combustion  rate. 


Tho  saaimum  amount  of  reaction  which  might  occur  between 
impinging  hydrazine/nitrogentetroxide  streams  was  estimated  using 
a  conservative  mixing  model  and  the  heat  release  data  of  Feiler  and 
Somogyi  (Ref.  13).  In  conjunction  with  this  nixing  model,  criteria 
was  established  to  predict  the  pressure  potential  tending  to  . 
separate  a  nixed  hypergelic  emulsion.  The  percentage  of  liquid 
phase  reaction  and  tho  subsequent  separation  of  tho  unlike  hyper* 
golic  streams  was  thus  predicted  for  tho  hydrasine/NTO  system. 

Tho  physical  model  and  tho  available  data  also  lead  to  the  postu* 
latios  of  an  ignition  index  numbor  suitable  for  rating  the  start* 
up  and  operational  characteristics  of  hyporgolic  liquid  phase 
contacting  devices, 

1.  Liquid  Phase  Reaction. 

a.  Heat  Release  Ratos. 

Heat  release  rates  of  hypergelic  liquid  systems  have 
been  measured  by  Feiler  and  Somogyi  (Ref.ig).  Their  experimental 
procedure  consisted  of  a  tangential  mixing  injector  and  a  nixing 
cup  of  length  L.  The  reaction  was  quenched  as  soon  as  it  left 
tho  nixing  length  L.  Tho  heat  release  results  of  these  experi* 
ments  aro  shown  in  Figure  73. 

The  shape  of  these  curves  aro  characteristic  of  the  degree 
of  nixing,  F(M),  occurring  within  the  mixing  length  L.  For 
complete  mixing  the  data  determine  the  maximum  rate  of  the 
liquid  phase  heat  release,  Feiler  and  Somogyi  explain  tho 
shape  of  these  curves  in  the  following  way  (Ref.,  ig) . 

"The  increase  in  heat-release  rate  with  injection  velocity  at  the 
lower  velocities  (Fig..73)  may  thus  bo  attributed  to  an  increase 
in  intorfacial  area  or  in  F(M),  The  rapid  rise  in  heat-release 
rate  observed  at  injection  velocities  of  120  to  140  feet  per 
second  is  attributed  to  an  increase  in  temperature.  In  this 
velocity  range  it  appears  that  heat  generated  in  the  inter¬ 
facial  area  can  no  longer  be  dissipated  to  the  surrounding 
fluid  at  a  sufficient  rate,  and  therefore  the  temperature  of 
tho  reacting  material  increases,  which  results  in  an  exponential 
increase  in  hoat-roloaso  rate.  At  velocities  greater  than  about 
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Injoction  velocity,  V^,  ft/sec. 

FIGURE  73  HSAT-RSLEASE  RATE  AS  FUNCTION  OF  INJECTIO 
VELOCITY(Rof,  18). 


140  feet  per  second  the  he&t>release  rates  attain  constant 
values  or  perhaps  decrease  slightly*  which  indicates  that  tempera¬ 
ture  has  attained  a  constant  value.  This  temperature  is  probably 
a  characteristic  temperature  of  the  propellant  system  such  as 
the  flame  temperature.  It  is  possible  that  gas  evolution  rates 
accompanying  the  reaction  are  sufficient  at  the  highest  injec¬ 
tion  velocities  to  inhibit  the  mixing  process.  The  resulting 
decrease  in  F(M)  (interfacial  area)  would  thus  account  for  a 
decrease  in  the  heat-release  rate." 

The  maximum  heat  release  rates*  with  complete  mixing,  were 
thus  determined  to  be  83,000  and  48,000  (kcal)/(sec)  (mole  HNO.) 
for  hydrazine  and  UDMH,  respectively,  (Figure  73). 

b.  Interfacial  Resistance  to  Reaction. 

Experiments  (Ref.  18)  determined  that  the  maximum  liquid  phase 
heat  release  rate  was  independent  of  oxidant  to  fuel  weight  ratio 
(i.e.,  concentration),  as  shown  in  Figure  74.  It  was  postulated 
rather  that  the  maximum  heat  release  rate  was  dependent  upon  the 
stoichiometric  ratio  of  reactants  for  the  following  reasons:  1) 
the  rate  of  heat  release  appears  to  be  independent  of  concentration, 

2)  four  times  as  many  molecules  are  required  to  oxidize  a  molecule 
of  UDMH  as  are  required  to  oxidize  a  molecule  of  hydrazine,  and 

3)  the  heat  of  reaction  per  mole  of  acid  for  fuel  is  approxi¬ 
mately  the  sane  for  either  fuel  (UDMH-158,  H-147,  kcal/mole) . 

This  treatment  of  the  data  may  be  justified  by  considering 
an  interfacial  resistance  mechanism.  The  following  two  observa¬ 
tions  lead  to  the  conclusion  that  reaction  kinetics  are  not  the 
controlling  mechanism  in  liquid  phase  fixing: 

(1)  The  observed  rate  of  reaction  is  independent  of 
concentration  (Figure  74),  and 

(2)  Ihe  half-life  of  a  mixture  of  hydrazine/NTO  is  on 
the  order  jf  10”'  seconds,  according  to  gas  phase  kinetics  of 
Glass  and  Sawyer,  (Ref.  19).  This  kinetic  rate  is  more  than 
four  orders  of  magnitude  faster  than  the  liquid  phase  mixing 
reaction  rate  data  of  Figure  75.  Thus  a  mechanism  other  than 
kinetic  reaction  rate  must  be  controlling  the  rate  of  liquid 
phate  reaction.  The  mechanism  postulated  here  is  one  of  inter- 
facial  resistance. 

It  has  been  previously  pointed  out  that  hydrazine/NTO 
appear  to  be  immiscible  (Ref.  20).  If  the  surface  reaction 
tine  is  faster  than  the  characteristic  mixing  tine  and,  if 
there  is  sufficient  pre'jsure  potential  in  the  products  to 
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FIGURE  74^  EXPERIMENTAL 
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Oxidant-fuel  weight  ratio 
Injection  velocity  170  ft/sec 
1.50  ^  Nitric  acid 
2.25  and  hydrazine 

Nitric  acid  and 
1.2s '  UDMH 

2.87*  _  _ 


. . . .  ^ 
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iEAT-RELEASE  CURVES  POR  VARIOUS 

>*(Rof,  18) 
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aaintain  separation,  then  this  icaiseibilit^  @ight  be  due  to 
the  fact  that  a  drop  of  one  eonponent,  having  its  Qm  surface, 
uill  tend  to  t  intain  this  surface  uhen  contacted  with  the  other 
eonponent.  Under  those  conditions  an  interfacial  resistance  to 
reaction  will  be  established  and  will  produce  an  effect  analogous 
to  ieniscibility.  Conparison  of  reaction  tine  to  raining  ti^e 
indicates  this  interfacial  resistance  nechanisn,  while  the  condition 
of  pressure  potential  is  treated  below. 

Considering  this  interfacial  resistance  to  reaction,  the 
rate  of  reaction  will  be  controlled  by  the  stoichiometric  ratio. 

This  is  duo  to  the  fact  that  the  two  liquids,  facing  each  other 
across  this  resistance,  are  both  pure.  Thus  reaction  is  dependent 
on  tho  stoichiouotric  ratio  of  each  fuel  systen  and  the  saxinus 
rate  with  perfect  nixing  will  give  a  conservative  estimate  of 
heat  release  for  any  hypergolic  fuel  systen,  based  upon  its 
stoichionotric  ratio  and  the  assunption  that  its  reaction  rate 
is  slower  than  its  nixing  rate.  Figure  75  presents  an  estimate 
of  tho  naxinun  heat  release  rate  produced  upon  conplete  nixing 
of  hypergolic  propellants.  This  naxinun  rate  is  valid  as  long  as 
the  propellants  are  confined  together  and  are  perfectly  nixed. 

c.  Characteristic  Confined  Mixing  Length. 

Tho  physical  aodel  of  Figure  76  was  postulated  in 
order  to  estinate  the  length  of  tine  during  which  inpinging 
streams  are  hold  together.  An  elenent  of  fluid  is  within  the 
inpingonent  area  for  an  average  length  of  tine, 

vD  , 

T  .  — 2.  1- 

^nix  8  v^  • 

where:  0^  is  the  jot  dianeter,  and  v^  is  the  forward  velocity 

of  the  jot  as  shown  in  Figure  76.  While  the  elenent  of  liquid 
is  within  this  area  of  inpingonent  it  is  confined  in  the  sane  way  as 
if  it  wore  contained  within  a  nixing  cup  sinilar  to  the  experiment 
of  Reference  18.  Although  perfect  nixing  nay  not  occur  for  in-  ^ 
pinging  stroans,  as  in  the  case  of  the  nixing  cup,  the  assumption 
of  perfect  nixing  sots  a  naxinun  limit  upon  the  heat  generation 
rate.  Thus  tho  average  tiao  within  tho  area  of  inpingonent 
gives  a  conservative  estinate  of  the  amount  of  forced  nixing 
occurring  between  two  unlike  hypergolic  streams.  Vihatever 
happens- after  tho  olooont  leaves  this  area  of  influence  will  be 
independent  upon  a  force  balinco  between  the  elenent *s  nonentun 
and  the  pressure  potential  of  reactant  gases.  This  problem  of 
stroan  separation  is  discussed  in  the  next  section  of  this 
report,  in  this  section  tho  purpose  is  to  calculate  the  naxinun 
aQOunt  of  propellant  reacted  during  liquid  phase  nixing. 
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f^oid  Retie  of  Reactants 


FIGURE  75.  MAXIMUM  HYFSRGOUC  HEAT^EELEASE  RATE, 
FOR  HYDRAZINE  TYPE  FUELS,  DATA  FROM 
FEILER  AND  SO^OGY!  {Rof.  18). 


10/ 


d.  Nunerical  Exanpie. 

The  following  exiople  is  based  upon  the  Bodel  of 
Figure  4  and  the  naxinun  heat  release  rate  data  of  Figure  7S. 
Consider  the  unlike  hypergolic  inpingenent  of  hydrazine/nitrogen- 
tetroxide.  Typical  paraneter  values  are:  jet  dianeter  equal  to 
0.050  inch  dianeter;  jet  velocity  equal  to  720  inches  per  second; 
and  inpingenent  angle  equal  to  90*.  The  characteristic  tine 
of  nixing  is 

1 

^nix  *  “s'”  v^  * 


The  stoichionetric  ratio  of  reactants  for  this  systen  is 
O.S.  Thus  fron  Figure  75  the  estinated  naxinun  heat  release 
rate,  with  perfect  nixing  is  62,000  cal/(sec)  (nole  acid). 

y  QS  A 

The  i^eat  generated  by  liquid  phase  reaction  is  given  as 


‘nix 
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ttlTT 


14)_(SxlO“'b 


4x10' 


Qn 
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nax 


sec. 


nix 


Thus 


Qliq  ■  24  cal/nole  of  oxidizer 


Whereas  the  heat  of  reaction  which  will  finally  bo  released  for 
the  hydrazine/nitrogen  tetroxide  systen  is  ^  >  248,000  cal/ 

nole  oxidizer.  Thus  the  anount  of  heat  roleaseo^ouring  liquid 
phase  nixing  of  a.typical  unlike  inpinging  doublet  is  estinated 
at  0.01%  or  1x10“^  fraction  of  the  total  anount  of  heat  available. 

A  sinilar  inpinging  MMH/NTO  doublet  will  have  the  sane  value  for  T  . 
while  its  reactant  nole  ratio  of  0.8  gives  a  naxinun  rate  of  heat 
release  of  83,000  cal/sec.  Thus  the  estinated  liquid  phase  heat 
release  is  33  cal/nole  of  oxidizer  and  the  percentage  reacted 
based  upon  a  total  heat  release  of  218  kcal/nole  NTO  is  0.015%. 


2.  Liquid  Phase  Separation  Potential. 

When  chenical  reactien  occurs  there  is  a  tenperature  poten> 
tial.  developed  at  the  reacting  interface  which  is  characteristic  of 

a.  The  nagnitude  df  the  heat  of  reaction. 

b.  The  enthalpy  capacity  of  the  innodiate  surroundings. 

c.  The  thernal  transport  properties  of  the  innediato 
surroundings. 
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If  conditions  of  eonfinonent  and  transport  properti@s  are 
dofinod  in  tho  icnedi^to  surroundings »  a  physical  nodel  nay  be 
postulated  so  that  absolute  nuebors  for  the  transient  tomperao 
turo  rise  nay  be  calculated.  However,  if  the  innediate  surround¬ 
ings  are  sooouhat  undefined,  an  index  nusber  nay  still  be 
dovolopod  by  caking  use  of  known  characteristics  of  the  systera. 
For  tho  case  of  reaction  between  icpinging  hypsrgolic  streams, 
tho  ieoodiat®  surroundings  are  not  well  defined. 


In  this  section  an  ostication  of  icpinging  stream  separation 
was  developed  for  two  cases. 

First,  for  tho  adiabatic  case  shich  neglects  thermal 
transport  in  the  surrounding  liquid,  and 


Seconds  for  the  unknown  thorcal  case  which  requires 
definition  of  an  index  nucbor,  considering  that  the 
thermal  transport  will  be  characterized  by  this  number 
and  thus  that  the  tendency  of  streaas  to  separate  may 
be  rated  according  to  this  index. 

a.  iiypergoiic  Stream  Impingement. 


The  two  stroams  of  Figure  4  impinge  at  2$  degrees  on 
a  comoon  plane.  The  common  inpingoment  area  is 

A  «  area  of  ellipse. 


Tho  normal  force  of  impingement  is  determined  from  a  momentum 
balance  on  tho  stream  before  and  after  impingement: 

-  wv 


n 


Sc 


sin  ^ 


This  force  results  in  a  pressure  at  A  such  that, 

F„  ■  AP  X  A.  HattJ  Rup©  (Ref.  ?1)  has  treated  the  dynami 
characteristics  of  free  liquid  jots  in  detail.  However,  detail  has 
boon  omitted  from  this  ordor-of-®agnitude  analysis. 

substituting: 


8c 


For  impinging  reactive  liquids  the  teaperatwr®  rise  of 
product  vapor  is  given  by  a  heat  balance 


AT  « 


where 


•  heat  of  reaction 
AHy  •  heat  of  vaporization 


aH|^  >  heat  lost  to  the  surroundings 


This  temperature  potential  may  he  related  to  the  pressure 
potential  necessary  to  overcome  the  momentum  of  impingement 
by  the  Clapeyron  equation 


ap 


/  A^^ 

'*  gl 


AT 


'AH„P>  . 

V  -at. 

vnRT2/ 


where  the  coefficient  (AH^/T  Vg)  is  the  average  slope  of  the 
vapor  pressure  curve  in  the  range  of  interest* 

Substituting  the  temperature  potential  for  pressure  results 
in  the  criteria  that: 


If, 


AHj^  •  AHy  "Allj^ 
Cp 


nRT 


4W  y  (sin 

’oj  Ic 


the  streams  will  be  blown  apart  by  reaction  and  thus  cannot 
be  mixed  by  impact  momentum* 


b.  Adiabatic  Case* 

If  it  is  assumed  that  no  heat  is  lost  to  the  surround! 
liquid,  (i.e,,  the  reaction  time  is  small  compared  to  the  time 
required  for  thermal  transport  »  then  equal  to  zero* 

Thus  if 


'  Allp  -  AHy  nRT^  4Wv(sin 
Cp  ^  AHyP  itDj 

the  streams  will  be  blown  apart* 

An  example  is  made  of  the  hydrazine/NTO  system  considered 
previously*  The  numerical  values  are 

AHj^  a  AHy  82*6  kcal/mole  of  reactant 

Cp  ■  8*6  X  10“^  cal/(molc  of  products)  (®K) 


with  7  moles  of  products  per  3  moles  of  reactant. 


Thus  the  adiabatic  toaperature  potential  (evaluated  at  200*C} 
available  in  this  systeo  is 


■  7  rrstijrs  •  «»*»'>* 

Whereas  the  tenperaturo  potential  to  be  ovetcone  in  order  to 
aoparato  the  streans  of  Figure  76  with  nuoerical  values  of 
oxaople  in  paragraph  lV>l.>d.,  evaluated  at  200*C  and  300 
pti,  is 

Teppcrature-Pressure  Conversion 

■  3.9  X  10^^  ‘K/atn 

Monontuo  Pressure 

if  .  =  1.65  aw 

Thus 

nM  .  6.42x10*  'K 

"v  •Bj  «c 

4  2 

Since  4,1x10  >  6.4x10  ,  the  hydrazine/NTO  streans  nay  be 

predicted  to  separate  after  initial  contact  and  surface 
reaction. 


c.  Separation  Index  Nunber. 

Since  the  stroan  inpingonent  is  not  truly  adiabatic 
a  separation  index  based  upon  the  Clapyron  potential  conversion 
can  be  used  to  evaluate  hyporgclic  strean  separation. 

This  index  is 


Is  • 


(AH^-  AH^) 

5 


nRT 

ah;^ 


4Wv(sin  6)^ 


Sc 


s©  that  tho  larger  the  value  of  I^  tho  greater  the  tendency 
of  the  stroaos  to  aoparato. 


3. 


Characteristic  Ignition  lndex« 


The  data  reported  by  Feller  §  Somagyi  (Ref.  is}  t<'hich  is 
presented  in  Figure  7:^,  can  be  used  to  develop  an  index  to 
relate  ignition  characteristics  of  different  injector  devices. 

The  words  which  Feiler  uses  to  describe  the  rapid  rise  in  heat 
release  rate  observed  in  Figure  73  can  be  interpreted  as  defining 
the  beginning  of  ignition  (Ref.  18): 

"The  rapid  rise  in  heat>release  rate  observed  at  injection 
velocities  of  120  to  140  feet  per  second  is  attributed  to 
an  increase  in  temperature*  In  this  velocity  range  it 
appears  that  heat  generated  in  the  interfacial  area  can  no 
longer  be  dissipated  to  the  surrounding  fluid  at  a  sufficient 
rate,  and  therefore  the  temperature  of  the  reacting  material 
increases,  which  results  in  an  exponential  increase  in 
heat- release  rate." 

The  onset  of  ignition  in  the  terms  defined  by  Feiler  repre> 
sents  a  microscopic  view  of  the  heat  release  per  unit  of  heat 
available.  Macroscopical ly  this  heat  released  determines  the 
average  rise  in  temperature  of  the  liquid  (i*e.,  after  thermal 
transport  has  occurred.)  The  heat  is  released  on  a  microscopic 
scale  at  an  immisible  interface,  causing  local  heating.  When 
sufficient  contact  time  is  allowed,  heat  is  available  to  raise 
the  average  fluid  cemperature  to  an  ignition  level.  The  per¬ 
centage  of  liquid-phase  heat  released  is  thus  a  guide  to  ignition. 
It  should  prove  that  when  this  number  exceeds  a  certain  value 
smooth  ignition  will  occur.  The  percentage  of  heat  released 
with  impingement  should  thus  prove  to  be  a  useful  ignition 
indexing  system  for  rating  the  smoothness  of  start-up  and 
operation  of  injector  devices.  This  index  number  is  given  by 
the  formula: 

n  T  . 

T  .  "max  mix 

where  R^j^jj  (Figure  7S  )  and  ar®  defined  s  in  Section  IV-1. 
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Conclusions 


(a)  The  percentage  liquid  phase  heat  release  of  typical 
iiapinging  hypergolic  streams  (MMil,  Hydrazine/NTO)  was  estimated 
at  0,01%.  If  this  amount  of  heat  went  to  heat  the  remaining 
liquid  it  would  lead  to  a  temperature  rise  of  approximately 
4*K,  The  percentage  of  liquid  reacted  upon  impingement  appears 
to  be  negligible,  while  the  amount  of  heat  released  may  cause 

a  measurable  liquid  phase  temperature  rise.  However,  the  cal¬ 
culated  4*K  temperature  rise  is  not  significant  enough  to  warrant 
its  inclusion  in  state-of-the-art  rocket  chamber  calculations, 

(b)  The  data  used  to  estimate  the  percentage  liquid  phase 
heat  release  also  lead  to  the  postulation  of  an  interfacial  resis¬ 
tance  to  mass  transfer  as  the  reaction  rate  controlling  mechanism. 
It  was  postulated  that  this  interfacial  resistance  would  pre¬ 
dominate  when  reaction  kinetics  were  faster  than  mixing  times 

and  when  there  was  sufficient  pressure  potential  generated  from 
the  reaction.  These  two  conditions  were  satisfied  for  the 
cases  of  Ilydrazine/NTO  and  MMH/NTO. 

(c)  Therefor.',  based  upon  conclusions  (a)  and  (b),  typical 
MMII/NTO  or  Ilydrazine/NTO  streams  will  separate  without  significant 
heat  release  under  conditions  of  simple  impingement. 

(d)  Feiler  and  Somogyi  explain  their  heat  release  curves 
in  terms  of  a  tendency  to  local  ignition  at  the  interfacial  sur¬ 
face.  based  upon  their  physical  description  and  their  measurement 
of  maximum  heat  release  rates,  an  index 

for  rating  ignition  and  start-up  devices  was  developed. 
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V.  ATOMIZATION 


Combustion  calculations  based  upon  droplet  vaporization 
models  require  that  the  distribution  of  drop  sizes  be  known# 

Many  experiments  have  been  performed  and  empirical  distribution 
functions  have  been  determined.  However,  there  is  little  current 
agreement  on  which  distribution  function  is  most  satisfactory. 
Since  all  the  empirical  distribution  functions  have  particular 
shortcomings,  a  theoretical  distribution  function  has  been 
derived  here.  Among  other  advantages,  this  theoretical  dis¬ 
tribution  allows  relation  of  the  droplet  size  to  the  physical 
injection  process.  Derivation  of  this  distribution  involves 
determining  the  most  probable  way  in  which  a  given  amount  of 
surface  energy  would  be  distributed  among  a  group  of  drops 
having  a  given  total  mass.  The  derivation  is  presented 
here;  and  it  is  shown  that  the  theoretical  distribution  fits 
experimental  data  and  can  be  related  to  atomization  efficiency. 

Of  the  several  better  known  distribution  functions 
commonly  used  to  describe  atomized  drop  sizes  none  have  been 
related  by  theory  to  parameters  of  the  atomization  process. 

The  Nukiyama-Tanasawa,  (Kef.  22)  Rosin-Rams ler,  (Kef,  33), 
Logarithmic-normal,  (Ref.  24  )and  Upper-Limit  distributions 
(Ref.  25  )  ore  primarily  useful  for  fitting  curves  to  observed 
data.  The  first  three  contain  two  distribution  parameters, 
one  representing  some  mean  diameter  and  the  other  the  dispersion. 
The  values  of  these  parameters  are  adjusted  to  fit  the  data  in 
question.  The  last  function,  the  Upper-Limit,  is  a  modification 
of  the  Logarithmic-normal  distribution  to  give  the  distribution 
an  upper  limit  to  confori><  better  with  some  observed  data.  The 
former  three  distributions  allow  for  infinite  drop  diameters  as 
a  possibility.  The  Upper- limit  distribution  includes  a  third 
parameter,  which,  however,  may  be  related  to  the  conditions 
of  the  atomization  process.  At  best,  then,  with  any  of  these 
distribution  functions  the  mean  drop  size  and  dispersion  can¬ 
not  be  characterized  except  through  observation.  However, 
there  is  a  distribution  that  might  be  related  to  the  atomization 
process.  That  would  be  the  most  probable  thermodynamic  dis¬ 
tribution,  Attempts  to  apply  a  thermodynamic  distribution 
CO  atomizew  dispersions  has  n  .t  been  reported 
and  should  be  examined. 


1 .  Theoretical  Development 


Consider  an  atoni^er  or  nozzle  which  in  a  given  increment 
of  tine  discharges  total  number  of  drops,  having  associated 
with  tho.9  a  quantity  9f  energy  Ef,  where 

Nt  *  f 

an^ 

:  E?  •  j 

where  ■  number  of  drops  possessing  a  particular  energy 
level/drop« 

^i  *  energy  level  of  a  single  drop  belonging  to  the 

i^*'  group. 

If  all  energy  states  have  equal  probability  of  occurring,  and 
the  drops  are  distinguishable  with  no  degeneracy,  with  the  two 
constraints  of  and  Ef  being  fixed  the  most  probable  thermo* 
dynamic  distribution  is  the  Maxwoll-Boltzman  distribution 

Nj  - 

whore  A  and  0  are  Lagrangian  multipliers. 

Then 

Ej  •  I  ci  o“^“®®i 


Noh  if  the  energy  of  each  drop  is  related  to  the  size 
of  the  drop,  the  drop  size  distribution  can  bo  determined.  For 
example  suppose  that  the  distributed  energy  is  that  of  surface 
energy,  then, 

o  oftD^ 

whore  a  is  the  energy  per  unit  surface  area  of  the  drop. 


For  largo  numbers  of  drops  we  can  write. 


N'l' 


II 


and 


and  the  probability  density  function  associated  with  F(v) 
is 


Furthernore  the  total  volume  Vj  of  the  Nf  drops  can  be 
shown  to  be 

Vt  «  ^  »2fl  Ny  o|o 
which  since  D3Q  is  defined  by 

Vt  •  ^  Nj  sDjJ 


leads  to  a  relationship  between  D20  Ojq. 
3 


From  this  the  Sautor-mean  diameter  can  be  found  to  be 


D32 


D202 


i  D 

v¥  *^20 


and  the  ratio  of  total  energy  to  total  mass 


^  60 5  _o 

^’'^T  "^^32  v¥  0^20 

The  last  expression  is  a  most  significant  one  in  that 
it  rebates  tho  ratio  of  total  energy  per  mass  to  the  average 
surface  diameter  of  the  atomized  spray  through  the  physical 
properties  of  the  fluid  system.  Thus  if  the  energy  per  mass 
could  bo  dotormined  from  the  physical  circumstances  for  a 
particular  atomization  tho  drop  size  distribution  would  be 
completely  described. 

Up  until  now  no  limitations  have  been  placed  on  tho 
possible  dianotors  of  tho  drops.  However,  another  distri¬ 
bution  can  be  dovolopod  by  an  identical  approach  which  limit 
tho  possible  diamotors  to  a  range  between  a  lower  and  upper 
dioQotor  Dj,  and  Dy,  The  resulting  equations  arc, 


F(N) 


erf{/P)  -  erf(,>T>j^) 


and 


erfC/Pu)  -  erf(/Pi^) 
where  P  «  ^ 


"I 


where  y(V(X)  is  the  incenplete  gamma  function  defined  as 
y(v»x)  » 


f  t'’-*  e-‘  dt 


'W 

Note  that  A  erf(x)  ■  y(4»x^'|  so  that  F(N)  can  b®  written 


finally  as  Cumulative  nuriber 'distribution 


F(N) 


.  ^  I,’’’) . :ii 


and  the  Probability  density  becomes 

f(N)  -  1 — ^224^ 


/2 


.-p 


The  remaining  final  equations  corresponding  to  those  developed 
for  no  maximum  or  minimum  diameter  are 


Cumulative  volume  distribution 


F(V)  « 

HPu)  -  HPl) 

Volume  probability  density 


f(V)  ■  2^2 


c  ’D20/ 


'1'  (Py)  - 
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^oittionsKip  between  D20  D3Q 


“ioV  .  »(Pu)  •  »(Pi,) 

®20  /  Y  P  \  P  ' 

If  T-'u; 

Sauter  oeen  dianoter 

D32  ■  020  )•*“’■•) 

t(y**’u)l  v(f.'’i) 

When  Pi,  G  &nd  Py  *  these  equations  reduce  to  the  simpler 
ease  given  before* 

2*  Theoretical  Distribution  Function  Compared  to  Experimental 
Bata.  '  '  '' . ''  ''  ' 


At  the  present  time  only  a  preliminary  evaluation  of  the 
usefulness  of  the  derived  distribution  has  been  made.  Experi¬ 
mentally  measured  distributions  by  several  investigators  have 
boon  compared  to  the  distribution  fuvictions  above  as  shown  in 
Figures. V- 1,2, 3Th93o  comparisons  pertain  only  to  the  ability 
of  the  functions  to  account  for  the  observed  dispersion  of 
drop  diameters  purely  as  a  function  of  the  observed  mean 
drop  dianoter,  in  this  case  the  surface  average  diameter. 

In  other  words  the  comparisons  are  merely  tests  of  the  distri¬ 
bution  function  to  '*fit"  the  data* 

The  results  of  those  first  conpariso.'is  are  very  encouraging. 
In  some  cases  the  theoretical  distribution  does  an  excellent 
job  of  predicting  the  number  of  drops  even  at  the  extremes 
of  the  diameter  ranges.  Furthermore,  in  situations  whore, 
only  some  moan  diameter  wore  known  and  nothing  known  about 
the  dispersioUf  the  theoretical  distribution  might  provide 
a  reasonable  ostinato  of  the  distribution  for  many  engineering 
problems.  (Note:  a  moan  diameter  based  upon  sny  definition 
is  suitable  since  for  a  given  distribution  function  all  moan 
diamotors  can  bo  related).  Nhore  this  is  possible  the  theore¬ 
tical  distribution  has  a  distinct  advantage  over  the  empirical 
distributions  which  havo  been  proposed,  since  various  investi¬ 
gators  have  dotoroinod  empirical  models  for  predicting  mean 
diaootors  for  some  tioszl«  types,  but  in  many  cases  have  failed 
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to  relate  the  drop  size  dispersion  to  any  useful  physical 
perasetersp  and  therefore  cannot  deteraine  the  distribution 
where  needed* 

A  second  significant  feature  of  the  theoretical  distri*' 
bution  is  that  it  opens  an  avenue  of  investigation  for  pre» 
dieting  drop  sizes  heretofore  only  briefly  explored*  Accord* 
ing  to  the  developcont  of  the  distribution  function  the  nean 
surface  diacoter  jhould  be  related  to  the  specific  surface 
energy  of  the  spray*  This  ieplies  that  exaaination  of  the 
energy  exchanges  involved  in  the  atonizati^n  process  of  soae 
nozzles  sight  be  used  to  predict  drop  aize  distributions* 
thus  a  thereodynasic  approach  rather  than  a  fluid  stchanical 
one  say  prove  fruitful.  Certainly  for  soso  nozzles^  forshla* 
tion  of  a  fluid  nechanical  nodel  is  alsost  hopelessly  eosplex, 
when  for  the  sane  nozzle  a  thersodynasic  exasination  is  feasible* 

In  the  developsont  of  the  distribution  function 
prosonted  here  only  the  surface  energy  of  the  drops 
was  considered*  Actually  other  forss  of  energy  such  as 
kinetic  and  vibrational  energies  sight  bo  distributed  as  a 
function  of  drop  size  and  would  affect  the  thersodynasic 
distribution*  Soso  consideration  of  these  factors  should  be 
included  in  any  future  investigations  of  the  relation  between 
drop  size  distribution  and  the  distributed  energy. 

3*  Correia- ion  of  Atosization  to  Process  Parasetors* 

Considering  only  the  surface  energy  it  has  boon  shown 
that  the  sprays  fros  sose  nozzles  can  be  '*fit''  very  well 
by  the  theoretical  distribution.  Two  questions  for  issediate 
investigation  would  appear  to  be: 

(a)  Can  the  deviations  of  soso  distributions  frosa  tho 
theoretical  distribution  be  related  to  nozzle  type*  experioental 
Bothod  used*  or  other  physical  paraaetors  of  the  syston? 

(b)  Can  the  surface  energy  of  tho  drops  forced  by  a 
nozzle  bo  calculated  froB  exaaination  of  tho  energy  exchanges 
in  the  atoaization  process  with  sufficient  accuracy  to  be 
useful  in  predicting  aoan  drop  diacoter?  If  so*  then  tho 
offieijney  of  an  atooizor  in  converting  kinetic  energy  into 
surface  energy  could  bo  correlated  directly  to  give  cean  dia* 
Botor  and  thus  droplet  distribution* 
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On  the  baais  of  vbat  has  been  aecoapliahed  so  far  it 
can  bO  tsaid  that  a  theoretical  distribution  of  drop  sises 
can.be  derived  nhich  is  related  to  the  physical  process  of  . 
atodiaatlen  and  the  properties  of  the  fluid.  These  parsaeters 
deteraine  the  efficiency  of  surface  energy  production  by 
parl^tular  atoaiter  devices*  Oreliainary  teats  of  this  dis¬ 
tribution  for  ^'fitting"  data  are  encouraging*  Ability  to 
predict  the  surface  energy  froa  the  energy  exchanges  in  the 
atoaisvtion  process  have  yet  to  be.  tested*  The  theoretical 
distribution  uhich  has. been  developed  here  should  serve  as  a 
valuable  liaiting  caso  lor  cosparison  and  help  in  correlation, 
in  addition  it  gives  physical  aeaning  to  distribution  fitting 
att^Qpts*. 
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bustion  problems .related  to  high  frequency  instability  in  liquid  rock* 
et  engines.  Using  the  steady-state  and  instability  computer  programs 
developed  under  this  study,  a  parametric  investigation  was  conducted. 

This  investigation  determined  the  influence  of  droplet  radius,  droplet 
distribution,  injection  velocity,  chamber  pressure,  and  mixture  ratio 
on  the  minimum  threshold  disturbance  required  to  trigger  combustion 
instability  in  a  Transtage  type  engine  configuration.  The  propellant 
combination  considered  was  monomethylhydrazihe/nitrogen  tetroxide.  Re.* 
suits  of  the  study  show  that  increases  in  injection  velocity  and  drop¬ 
let  distribution  increased  stability.  An  increase  in  chamber  pressure, 
based  on  constant  flow  rate  increased  stability  while  increases  in 
chamber  pressure,  at  a  constant  contraction  ratio,  decreased  stabiiit) . 
There  appears  to  be  a  droplet  size  for  minimum  stability,  with  change: 
in  either  direction  resulting  in  improved  stability.  Results  also  shot 
that  due  to  the  vapor  phase  reactions,  monpnethylhydrazine/nitrogen 
'jetroxide  vaporize  at  approximately  the  sane  rate.  Thus,  the  oxidizer 
or  fuel  could  be  made  to  control  the  combustion  process  by  slight 
changes  in  the  injector  and  engine  parameters.  For  the  engine  config-  I 
uration  studied  the  oxidizer  vaporized  slower  than  the  fuel.  A  pro-  I 
gram  review  summarizes  the  work  reported  in  the  Semiannual  (DSC  SN-68|S1) 
and  Special  report  (AFRPL-TR-65-2S4  written  under  this  program.  I 
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